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Abstract 

We investigate the multi-loop correlators and the multi-point functions for 
all of the scaling operators in unitary minimal conformal models coupled to 
two-dimensional gravity from the two-matrix model. We show that simple 
fusion rules for these scaling operators exist, and these are summarized in a 
compact form as fusion rules for loops. We clarify the role of the boundary 
operators and discuss its connection to how loops touch each other. We derive 
a general formula for the n-resolvent correlators, and point out the structure 
similar to the crossing symmetry of underlying conformal field theory. We 
discuss the connection of the boundary conditions of the loop correlators to 
the touching of loops for the case of the four-loop correlators. 
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1 Introduction 



Quantization of gravity is one of the most important issues in physics. The 
understanding of two-dimensional quantum gravity, which is the simplest quantum 
gravity, has experienced great progress through the study of matrix models J It was 
proposed [0 that the integral over the geometry of two-dimensional surface can be 
descretized as a sum over randomly triangulated surfaces and such regularized two- 
dimensional gravity can be realized by hermitian matrix models. Feynman diagrams 
of the matrix models correspond to the dynamically triangulated surfaces and the 
continuum limit of the models then describe the theory of two-dimensional gravity. 

Due to the double scaling limit [||] the sum of the contributions from all topologys 
of two-dimensional surface can be treated, and thereby the matrix models have been 
drawn much attention as a non-perturbative definition of non-critical string theories. 
Following the discovery of the double scaling limit, many important structures of 
the models have become clear; for example, the connection to KdV flow [ffl, the 
Virasoro and W constraints ^ |? 01- Field theory of non-critical strings ||, [10 
has been constructed based on the matrix models. 

The matrix models include infinite critical points, which are considered to repre- 
sent certain conformal matters coupled to two-dimensional gravity. The m-th criti- 
cal point of the one-matrix model corresponds to the (2m + 1, 2) minimal conformal 
model coupled to two-dimensional gravity. The general (p, q) minimal conformal 
model, where the central charge is c = 1 — 6(jl ~ g ) can be realized as the continuum 
limit of the discrete system where the degrees of freedom are points on the A(DE) 



Dynkin diagram ||12|| . Multi- matrix chain model has been introduced as a model 
which includes the critical points corresponding to the general (p, q) minimal mod- 
els coupled to gravity. In this model, q matrices interact as a chain. The two-matrix 
model O, [TJL |Tj|, which is the simplest multi-matrix chain model, however, turned 



out to include all (p, q) critical points, which was pointed out in JT3|, [L4| and shown 
explicitly in (|TjJ. We use the two-matrix model to investigate the unitary minimal 
model (m + l,m) coupled to two-dimensional gravity. 

The emergence of the infinite number of scaling operators is one of the most 
important properties of the matrix models. Before coupled to gravity, the minimal 
model has finite number of primary fields. Coupled to gravity, however, infinite 
number of scaling operators emerge. This phenomenon can be understood as fol- 
lows. In the Kac table we can divide the primary fields $ rjS into those which are 



3 See for example 0] for review. 

4 The corresponding structures in continuum framework have been shown also in 
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inside the the minimal conformal grid 1 < r < q — 1, 1 < s < p — 1 and those 
outside, which correspond to the null states. Before dressed by gravity, the fields 
outside the minimal conformal grid decouple |L7| from physical correlators. After 



gravitational dressing, they cease to decouple []lq , |19|j and become infinite number 
of scaling operators. The similar phenomenon has been shown in continuum frame- 
work. Through the examination of the BRST cohomology of the coupled system 
composed of Liouville theory, the ghosts and the minimal matter, infinite physical 



states were shown to exist [g0|, |21|| . These states have their counterparts in the 
matrix models as the scaling operators. Some of the scaling operators do not have 
their counterparts in the BRST cohomology, which we will discuss later. 

In ordinary (p, q) minimal conformal model the primary fields satisfy certain 
fusion rules [fL7H ; three-point function ($ ri Sl $ r2 S2 $ r3 S3 ) is non-vanishing only when 



1 + | r i — 7"2| < ^3 < mim{ri + r 2 — l,p}, T\ + r 2 + r% = odd 
1 + \s\ - s 2 \ < s 3 < mim{si + s 2 - 1, g}, Si + s 2 + s 3 = odd . (1.1) 

It is interesting to examine how the fusion rules change when the matter couples 
to gravity. In particular, we are interested in the fusion rules for the gravitational 
descendants (<jj, j = q + 1, q + 2, • - •), most of which correspond to the operators 
outside the minimal conformal grid. Before coupled to gravity the corresponding 
fusion rules do not exist. Such three-point functions were examined from the point 
of view of the generalized KdV flow in [B^] for lower dimensional scaling operators in 
the case of (m + 1, m) unitary matter. It was shown that the gravitational primaries 
o,j [j = 1, ■ ■ ■ , m — 1) satisfy fusion rules of BPZ type; for ji + j 2 + J3 < 2m — 1, 
((fj^jzVjz) is non- vanishing only when 

i + b'i-ial <h< ji+j'2-1 • (i-2) 



The fusion rules were also examined in continuum framework |Tj|. As for the grav- 
itational descendants, however, we think clear results have not been obtained. In 
this paper we would like to clarify the fusion rules for all of the scaling operators 
including the gravitational descendants in the case of unitary minimal model. This 



paper is based on E3L E4L 25, 26 



Macroscopic loop correlators ER E3, which are the amplitudes of the surfaces 



with boundaries (loops) of fixed lengths, are the fundamental amplitudes of the 
matrix models. Although these amplitudes are hard to treat in the continuum 
framework p2| , they are defined quite naturally in the matrix models. It was shown 



30| , |31| that these correlators have more information than those of local operators 
and that the latter correlators can be extracted from the former correlators explicitly 
in the case of c = 0,1/2,1. They argued there that macroscopic loops could be 
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replaced by a sum of local operators in a certain situation and thereby obtained the 
correlators of local operators from those of macroscopic loops. 

One of the purposes of this paper is to generalize the idea in [^0| to the cases of the 
general unitary minimal models and to clarify the fusion rules for macroscopic loops 
and all of the scaling operators. First, we derive the three- and n-loop correlators 
from the two- matrix model at the general unitary critical points pi, E5[, and then 



derive the explicit forms of the correlators of the scaling operators P5 |. The main 
conclusion is that the three-point correlators of all of the scaling operators satisfy 
certain simple fusion rules |2fJ and the fusion rules for all of the scaling operators 
are summarized in a compact form as the fusion rules for three-loop correlators . 

In matrix models, there are infinite subset of the scaling operators which do not 
have their counterparts in the BRST cohomology of Liouville theory. In the case 



of one-matrix model, Martinec, Moore and Seiberg |33[ argued that these operators 
are boundary operators, which correspond to the vertex operators of open string 
and couple to the boundaries of two-dimensional surface. They proved that one of 
them is in fact a boundary operator which measures the total loop length. We think, 
however, the roles of the rest of these operators were not clear. We also clarify the 
role of these operators and its connection to the touching of loops in the case of 
general unitary models p6| . 

We also determine completely the forms of the multi-resolvent correlators, which 
are the Laplace transform of the multi-loop correlators, and point out that the loop 
correlators have the structures similar to those of the crossing symmetry of the 
underlying conformal field theory |[25|| . In the cases of four- and five- loop, we discuss 



the connection of the boundary conditions of the loops to the touching of the loops 



26 



As another formulation of 2D gravity with matter system, models of string whose 
target spaces are the Dynkin diagrams have been investigated P4J . We also comment 
on the connection of our results to those from these models. 

The paper is organized as follows. Sect. 2 is devoted to the review of the matrix 
models and the macroscopic loops. We limit our discussion to the subjects that have 
direct connections to the subsequent sections. 

In sect. 3 we derive the three-loop correlators and extract the three-point func- 
tions for all of the scaling operators through expansion of loops in terms of the 
local operators. We then show that certain simple fusion rules exist for these local 
operators or loops. We also discuss the role of the boundary operators there. 

We derive the formulas for the multi-resolvent correlators in sect. 4, and give the 
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explicit forms of the four- and five-loop correlators. We point out that the structure 
corresponding to the crossing symmetry of the underlying conformal field theory 
exists in the multi-loop correlators. We also discuss the connection of the boundary 
condition of the loops to the touching of the loops. 

Sect. 6 is a summary. 
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2 Conformal field theory coupled to 
two-dimensional gravity (review) 



2.1 Matrix models and two-dimensional gravity 

Let us briefly review the matrix models and the connection to the Liouville theory, 
emphasizing on the notion of the scaling operators and that of the macroscopic 
loops. We limit our discussion to the subjects which have direct connections to the 
later sections. 



2.1.1 Matrix models and random triangulation 

Let us consider the model defined by the path integral with respect to an N x iV 
hermitian matrix $, 



d$e-5 ir * + ^ ir * , (2.1) 



where the measure is 



d$ = [] d$*i I] dRe^dlm^. . (2.2) 

i i<j 

The propagator is (^^^J = S^Sj, and is represented by the double lines in 
fig. ffl. (a). The arrows connect the upper matrix indices to the lower ones. The 
vertex in the action is represented by fig. [I] (b). Expanding the partition function in 
term of g, we find that the each Feynmam diagram represents a net on an orientable 
two-dimensional surface. Taking the dual of such a diagram, the vertices turn into 
triangles and the dual diagram represents a random triangulation of two-dimensional 
surface. Therefore, the model specified by eq. (|2.1| ) can be considered to represent 
a theory of random triangulation of 2D surfaces and is expected to be a theory of 
2D quantum gravity when we take continuum limit. 

Let us count the power of N associated to each diagram. Changing variables 
$ — > <&/y/~N, the action becomes iVTr (— | Tr $ 2 + g Tr $ 3 ) . From this form of 
action it is clear that each vertex contributes a factor of N, each propagator (edge) 
contributes a factor of AT -1 , and loop (face) contributes a factor of N due to the 
index summation associated. Each diagram has thus an overall factor 

N V-E+F = N X = N 2-2h ? (2 _ 3) 

where x an d h are the Euler character and the number of genera of the surface 
associated to the diagram respectively. From (|2.3|), the partition function can be 
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(a) 



(b) 



Figure 1: propagator and vertex 



expanded as 



Z{g) = Y,N 2 - 2h Z h {g) 



(2.4) 



where Zh(g) represents the contribution from the surfaces of genus h. In the large 
N limit, the contribution from the planar surfaces dominate. 



2.1.2 Continuum limit 

When expanded in the coupling g, for large order n, Zq behaves as 

Ua) ~ E^ 3 (-) ~ - gf- T - (2.5) 

n 9* 

and the expectation value of the number of vertices (triangles in the dual diagram) 
is given by 

(n) = ^-\nZ (g) — . (2.6) 

dg g-g* 

The partition function Z Q (g) thus becomes non- analytic and (n) diverges when g ap- 
proaches some critical value g*. Since (n) diverges as g — > g*, it is expected that the 
contribution from the continuum surface with finite area can be obtained by rescal- 
ing the area of the individual triangles to zero accordingly. The contribution from 
the continuum surface is considered to correspond the non-analytic part of eq. (|2.5| ) . 
Therefore, the behavior of the model (|2.1|) near the critical point is considered to 
represent two-dimensional quantum gravity. 
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2.1.3 Multi-critical points and multi-matrix models 



So far, we have considered the model( |2TT| ), which consist of one kind of vertex in the 
action. As a generalization of this model, let us consider the the model specified by 
the action 

N 

S = j 1W($) (2.7) 

where V(<$>) is some polynomial of matrix $. This model can represent a series 
of systems of consisting of matter and two-dimensional gravity. By tuning the 
couplings in the potential, various critical points are obtained. The m-th critical 
point corresponds to the (2m + 1,2) minimal conformal model coupled to two- 
dimensional gravity. 

As another generalization, let us consider the multi-matrix chain model, 

e z = ffl* a e- s , (2.8) 

v-l v-2 

S = J2 v «(® {a) ) - E c a $ (a) $ (Q+1) . (2.9) 

a=l a=l 

Here the different matrices <&( Q ) represent u — 1 different matter degrees of freedom 
that can exist at the vertices. Note the couplings c a in the action ( |2.9| ) couple the 
matrices along a line (chain). 

It was pointed out, however, in [y]] and shown explicitly infl5|] that it is 



sufficient to consider the two-matrix model in order to generate the most general 
critical points, which correspond to the c < 1 minimal conformal models. We will 
thus use the two-matrix model to examine the minimal models coupled to two- 
dimensional gravity in this article. 



2.2 Scaling operators in matrix models 

2.2.1 KdV flows and scaling operators 

Consider the two-matrix model with symmetric potential, 

e z = J dAdB e-£ ^ )+ u { b)-ab) ^ (2 10) 

where U is an arbitrary polynomial of order m. In this article, we limit our discussion 
to the two-matrix model with symmetric potential and to the critical points which 
correspond to the unitary minimal models. In the case of asymmetric potential, some 
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of the boundaries (loops) of the surface would have fractal dimensions different from 
the usual dimension of length. 



Integrating out the "angular" variables, we have 16 



e z = J dAjdAjA(A)A(A) e "f E,(^.W,)-a,a,) (2 _ n) 

Here A (A) is the Vandermonde determinant, A and A represent the eigenvalues of 
the matrices A and B respectively. 

We introduce the orthogonal polynomials \j) = £j(A) and (k\ = £fc(A) by the 
orthonormality relation 



(j\k) = J dfi &(A)&(A) = 5 jk , 

dfi = dXdl e-^ +u ^- x ~ x ) . (2.12) 

We define matrices A and P by 

A nm = (n\X\m) (2.13) 

Pnm = (n\ ^ \m) . (2.14) 



It is obvious from the definition ( |2.13 ) and Q2.14Q , that A and P obey the Heisenberg 
commutation relations 

[P,A] = 1. (2.15) 

Now the important fact is that the operators P and A have non-zero matrix elements 
Pij and Aij only if is sufficiently small. Since the bounds are independent of N, 
in the limit iV — > oo, P and A become differential operators (in x, the cosmological 
constant) of finite order. The continuum scaling limit of the two-matrix model is 
abstracted to the mathematical problem of finding solution to eq. ( j2.15|) . Let us 
consider the (m + 1, m) critical point which corresponds to the (m + 1, m) minimal 
model. After suitable renormalization, A is given by 

A = D m + u(x) m _ 2 D m - 2 + u(x) m ^D m - 3 + ■■■ + u(x) , D = d x . (2.16) 

(By a change of basis of the form A — ► f^ 1 (x)Af(x), the coefficient of D m ~ l may 
be always be set to zero.) and P is given by Q 

P = (L m+1 )+ . (2.17) 

Here L = A x l m is a pseudodifferential operator satisfying 

A = L m , L = D + a^- 1 + a 2 ZT 2 + ■ ■ • , (2.18) 
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and {L a ) + denotes the nonnegative (differential operator) part of L a . Substituting 
eqs. ( |2.160 and ( p. 171) into eq. Q2.15Q , we find that differential equations for u(x)i. 
These equations determine u(x)i up to m — 1 integration constants ij (i = 1, • • • , m — 
1). The dependence of L on the constants ti is given by the first m — 1 generalized 
KdV flows: 



{V 



(2.19) 



In general, the perturbation from the (m + l,m) point is represented by the inde- 
pendent flows in term of commuting operators {L % \i — 1, 2, • • • ^ {mod m)}\ 



d_ 

d_ 

dti 



A 
P 





' A 




~(L%, A 




, p 




"(L')_, P 



(2.20) 



The infinite number of directions of the perturbation correspond to flows along 
RG trajectories between various critical theories, identified with the (p, q) minimal 
model coupled to two-dimensional gravity. Since the perturbation is realized by 
adding an infinite number of relevant matter operators dressed by gravity to the 
original critical action, the correlation functions of the scaling operators are defined 
by the following relation: 



\(J h a j2 



a a ) 

Jn I 



_d d_ 



_d_ 

dU, 



\ogZ . 



(2.21) 



2.2.2 Correlators from KdV flow 



The correlation functions of the scaling operators <jj on the sphere were calculated 
in [[EJ for lower j from the point of view of the KdV flow. 

Let us define the following infinite number of commuting operators on the sphere: 



n 



3=0 



J 



with 



[Qn, Qk]=0 . 

The operators A and P on the sphere are given by 



(2.22) 



(2.23) 



.4 
P 



{Qm)+ — 
(Qm+l) + 



L m = D m _ _ muD m-2 + . 

2 



(2.24) 
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where u is the two-point function of dressed identity operator. Substituting eqs. 
fl23f) into eq. (]215D , we have 



u 
2 



x 



m + 1 



l/m 



(2.25) 



The correlation functions of Oj for lower j can be calculated from the KdV flow. 
For example, from eqs. ( [2.20] ), we obtain the following expression for the one-point 
function: 

^-u = -2(Res L n )' , (2.26) 

where Res L n is the coefficient of D~ x in L n and L n (n < 2m — 3) can be replaced 
with Q n due to the relation 



n 



Qn~L n = -c n D n - 2m + 0(D 
m 

The two-point functions can be calculated from 

d d r , B ... 



n—2m—2\ 



dt n dt k 



u = 2 (Res [LI, LI)) 



(2.27) 



(2.28) 



As for explicit results of the correlators, we mention these in sec. 3.5. 



2.2.3 Connection to continuum theory 

The m-th critical point of one-matrix model corresponds to the (2m + 1, 2) mini- 
mal model coupled to two-dimensional gravity. The scaling operators are naively 
expected to correspond to the following operators in the (2m + 1, 2) minimal model 
coupled to Liouville theory in the continuum framework: 



. «_> J e^fc^, . , j = 0, ■ • ■ , m - 2, (2.29) 

where atj = ^j(m — j), $ rjS are the primary fields of the corresponding conformal 
field theory. 



This correspondence fails however. In 30 , it was argued that the discrepancies 



were due to contact terms which arise when two operators are at coincident points. 
They showed explicitly the correct correspondence by the analytic redefinition of 
coupling constants tj 

t n = C\k + C%tij + ■■■ , (2-30) 

mainly for the case of one-matrix model. The original frame of operators <jj and 
couplings tj is referred to as the KdV frame and the new frame of operators dj and 
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tj is referred to as conformal field frame. The wave function of dj is proportional 
to the modified Bessel function K L (2JJit) so that it satisfies the (minisuperspace) 

a 

Wheeler-deWitt equation 



= , (2.31) 



which is a desirable property. 

The BRST cohomology of the coupled system of Liouville theory, ghosts and 
the (p, q) minimal matter was examined in |2U], |2T|]. It turned out that the BRST 



cohomology is spanned by infinite operators of the form 

Q. e <w °!± = P + q ~ 3 j > 1, ^ (modp), + (mod q) , (2.32) 



where <p is the Liouville field and 7 = (V25 — c — V 1 — c) / vl2. The operators Oj 
are made of ghosts, matter and derivatives of (p. On the other hand, the scaling 
operators cr, of matrix model at the (p, q) critical point scale like Liouville operators 
of the form 

Oj e a ^, ^ = P 3 j > 1, ^ (mod q) . (2.33) 

Apart from the discrepancy of the operators with j = (mod p), the two calculations 
are in remarkable agreement. It was argued in |33| that the scaling operators with 
j = (mod p) are boundary operators, which couple to the boundaries of two- 
dimensional surface and correspond to the vertex operators of open string. 



2.3 Macroscopic loops 

2.3.1 Macroscopic loops in two-matrix model 

In the two-matrix model, the operators 

Tri ni , TriT 2 (2.34) 

create holes with boundaries of lattice lengths n\ and ri2 respectively. The correlation 
functions of Tr A ni or Tr A ni are expected to become those of macroscopic loops in 
the limit ani — > with finite, when the unit lattice length a approaches zero. 

It is convenient to consider first the correlators of the resolvents 

W+( A )=Tr— , W-(p i ) = Tr^- 1 r, (2.35) 
Pi - A Pi - B 
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where pi is a parameter corresponding to the bare boundary cosmological constant 
of each loop. Due to the formal expansion 

00 Tr A n 

W + {Pi) = £ ' ( 2 - 36 ) 

n=0 Pi 

the resolvents include the contributions from loops of any length. The correlators 
become singular when pi approach some critical value p*. Since the contributions 
from loops of finite continuum length corresponds to those of infinite lattice length, 
continuum loop correlators are defined as the inverse-Laplace image of non-analytic 
part of the resolvent correlators with respect to Q = (pi — p*)/a. 



2.3.2 'Classical' solutions to Heisenberg relation 

In later sections we will use extensively the 'classical' solutions to the 'classical' 
Heisenberg relation. Let us explain these in this subsection p. 3.2 . 



Since we would like to examine the correlators on the sphere, we are interested 
only in the planar limit ( large N limit ). The Heisenberg relation ( 2.15Q turns into 



[^] = ^, (2-37) 

after rescaling P — > j^P- From eq. ( |2.37|) , we see that 4 plays the role of Planck 
constant. It is thus expected that the corresponding 'classical' functions would be 
much easier to handle than the operators A and P in the large N limit. 

At this point it is useful to change notation for the indices of the matrix elements: 

A k (n) = A n _ k}Tl , P k (n) = P n - k , n ■ (2.38) 

Here n represents the position of the matrix element on the diagonal, and k is its 
deviation from it. Then the action of the operators A and P on the orthogonal 
polynomial basis is described by 

m-l (m-l) 2 

A\n)= £ \n-k)A k (n) , P\n) = ]T \n — k) P k (n) . (2.39) 
fc=-i fc=-i 

The 'classical' functions are defined by 

m—l 

A(u,s) = £ e k "A k {n) 

k=-l 

{m-l) 2 

P(u,s) = £ e k "P k (n), (2.40) 

k=l 
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where s is the continuous variable 

s = —A (2.41) 

and us is its conjugate coordinate. 
The equation of motion 

^■ = (z|t/'(A)|i)-4, ( 2 - 42 ) 
which is obtained by doing an integration by parts, reads 

P{u>, s) = U' (A(u, s - . 1 _ A (-u + A^, s) . i , (2.43) 

when expressed in term of the classical functions. In the planar limit, eq. ( p.43| ) 

reads 

P(u, s) = U'(A(u, s)) - A(-w, s) , (2.44) 
and the Heisenberg commutation relation ([2.371) is replaced by the Poisson bracket 

{P{u , s) ,A(u,s)} = ----- = -l. (2.45) 

Note that in the large N limit the classical functions A and P depend on A only 
through s, which is easily seen from eqs. Q2.43|) and (|2.44j ). 

Let us find the solution to the classical Heisenberg relation eq. (|2.45|) near the 
(m+l,m) critical point, which corresponds to the (m + l,m) unitary minimal con- 
formal model. 

At the critical point, one expects the following singular behavior of A and P: 
A{z) - A, ~ (1 - z) m , P{z) -P.~(l- z) m+l . (2.46) 
From the scaling laws ( p.46|) , the solution to the Heisenberg relation is given by |15 



A(z, s) - A = 2r] m cosh m6 

P(z, s) - P* = 2r] m+1 cosh(m + 1)9 

s-A, = (m + l)r/ 2m . (2.47) 

Here P*, A* and A* denote the critical values of the corresponding quantities and 
the parametrization 

u = 2r]cosh6 (2.48) 

is used. We will use the classical functions ( |2.47|) extensively to calculate the loop 
correlators in later sections. 
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2.3.3 Loops in semi-classical Liouville theory 



When we discuss the loop correlators, it turns out to be very helpful to consider 
these correlators semi-classically in Liouville theory. Let us explain these |3(J briefly 
in this subsection |2.3.3| . In the continuum framework, two-dimensional gravity part 



of the coupled system can be described by Liouville theory based on the action, 
fi*33 = IT f d 2 ^g ab d a ^ 

+ ir([ d 2 ^% + / dsk<f) 

87T \JE V JdT, J 

+ ~J~2 I d 2 e^e>+ I dse^ 2 , (2.49) 

SlT'-f 2 Jt V 47C-f 2 JdT 

where g a b is a reference metric and g a b e 7 ^ is a physical metric, R and k are respec- 
tively the curvature and the extrinsic curvature of the boundary with respect to the 
reference metric g a b- We denote by p and p bulk and boundary cosmological con- 
stants respectively. Classically, Q = 2/7 where 7 is the Liouville coupling constant. 
Let us consider the correlation function 

(j[ e *Mzi)j = JVip e~ SL ]Je aMzi) . (2.50) 

We obtain the classical equation of motion: 

J_ A ^ -f e w + E a^ 2 \z -Zi) = 0. (2.51) 

47T 57 

Since the curvature of the physical metric is 

R[e^g] = -e-^A( W ) , (2.52) 



eq. (|2.51| ) describes a surface with constant negative curvature and the inserted 
operators e ai(p ^ play the role of the sources of curvature. Note that in the absence 
of a boundary a solution exists only when 

X = J2 a i + ?( 2/ >-2) ( 2 - 53 ) 

Li 

I 

is positive, where h is the number of handles. The nature of the surface and hence 
the nature of associated quantum states depend crucially on the sign of X. When 
there are boundaries, a classical solution always exists. Let us restrict our attention 
to the case with a single boundary and discuss whether the boundary can be replaced 
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by local sources of curvature. In this case the nature of the surface depends crucially 
on the sign of 

Y = X + ±Q = J2^-\qx- (2.54) 

Case 1 : Fixed /x, Y > 0. When the loop is shrunk to a point, there exists a 
classical solution with constant negative curvature. A small loop behaves like a 
local source of curvature Q/2. 

Case 2 : Fixed /x, Y < 0. When the loop is shrunk to a point, there is no classical 
solution with constant negative curvature. We can understand this case better if we 
constrain the area of the surface to be A. 

Case 2-1 : Fixed A » f 2 , Y < 0. The classical solution has positive constant 
curvature and the small-£ limit is smooth and the loop becomes as a puncture with 
curvature Q/2. 

Case 2-2 : Fixed A <C I 2 , Y < 0. The classical solution has negative constant 
curvature and the loop cannot be thought of as a local disturbance. 

So far we have discussed classically. Semi-classically, Y must be modified by 

(2.55) 

where a m i m is the curvature associated with the lowest dimension operator (9 m i n ( 
the dressed identity operator, in the case of the unitary minimal matter) because 
this is the maximum curvature that can be localized in a point in the quantum 
theory. Similar observations follow in semi-classical discussion. In case 1 and case 
2-1 the loop can be replaced by a sum of local operators and the contribution to 
the amplitude give rise to non-analytic terms in /x. In case 2-2 the loop cannot be 
replaced by a sum of local operators and the contribution to the amplitude give rise 
to analytic terms in /x; the loop length t plays the role of ultraviolet cutoff. 
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3 Three-loop correlators and fusion rules 



In this section we consider the loop correlators in the unitary minimal models 
(m + 1 , m) coupled to two dimensional gravity and the physical information we 
can extract from these. As shown in the case of one-matrix model [[!(]], the loop 
correlators are expected to have much more information than those of local opera- 
tors. We calculate the three-loop correlators in the systems stated above, from the 
two-matrix model with symmetric potential, at the (m + 1, m) critical points and 
show that simple fusion rules exist for the loop correlators and for all of the scaling 
operators. 

3.1 Formula for n-resolvent correlator 



Consider the connected part of the n-point correlators of the resolvents, which we 



introduced in sec. 2.3.1 



W+ {p .)= TY _}_^ ; w-( Pj ) = Tr—^ . (3.1) 
Pi - A pj- B 

First, let us show briefly the formula for the n-resolvent correlators we obtained 
in p5[i. The explicit derivation of the formula will be shown in section |4.1| later. At 



the (m + l,m) critical point, we obtained the following formula for the n-resolvent 
correlator: 

(Tl^^hni-ehy^^ ' (3 ' 2) 



Here we denote by ((• • •)) the connected part of the averaging with respect to the 
matrix integrations and is some function of Q and A, through z* and their 
derivatives with respect to the bare cosmological constant Aj. Note that we intro- 
duced independent cosmological constants Aj for each loop for the convenience of 
the calculation. We put Aj = A at the end of the calculation in eq. ( |3.2| ). The 
function z* of Q and Aj is parametrized as follows, 

z* = exp(2?7j cosh #j) , (3.3) 

where 

Pi-p^ = ad = A(z*; Aj) - A* = aMi cosh wBi , t) { = (aMj/2) 1/m , (3.4) 



5 The formula for the multi-loop amplitudes in the case of the general one-matrix model was 
derived in fl29f . 
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A; - A* 



(m + 1)77, 



2ni 



Mi 

- ,r l>i • \ ~2 



(3.5) 



m + 1 

where p* and A* represent the critical values of Pi and A(z; A) respectively. We 
denote by Q and /x, the renormalized boundary and bulk cosmological constants for 
the corresponding loop respectively. 



The origin of the parametrization eq. (|3.3| ) comes from the planar solution to the 
Heisenberg algebra ( |2.47|) . In fact, the function A(z, A) in eq. (|3~4l ) represents the 
solution at the (m + l,m) critical point. 

The function is easily written down for lower n. For n = 2, 3, we have 







dz* 



1=1 % j (A) 



Z; — Z 



dz. 



* 3 



e aa. n 

i=i ,-(?«) 



2,- — Z A 



(3.6) 
(3.7) 



For n = 4, 5, the correlators can be written compactly using graphs as introduced 
below: 



R^ 



E 



a 



OA,. 



14 




12 



+ E 



13 11 12 l 4 

o-#=»-»o 



(3.8) 



R^ 



E 



(9 \2 



9A n 



u o- 



nla 



61 5 



+ E 



9A n 




+ E 



14 h 12 13 15 



(3.9) 
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In these figures a double line linking circle i and circle j, a single line having an 
arrow from circle i to circle j and a solid circle % represent (z* — z j)~ 2 , ( z * — z j)^ 
and -q£- respectively. The summations are over all possible graphs that have the 
same topology specified. Each graph appears just for once in the summation. Note 
that the links to the external circles are not double lines but the single ones with 
arrows and that the internal circles are solid circles. 

For general n, the function is expressed in the same way. The rule is as 
follows. First, we consider all possible graphs which have n circles and n — 1 links in 
the same way as in the case n = 5. Second, if the internal solid circle i has l{ links 
in each graph, the graph is operated by IL(af~) * 2 ■ Then the summation over all 
graphs gives the expression for R( n \ 



3.2 Three-loop correlators 

3.2.1 Derivation of three-resolvent correlators 

As a example, let us calculate explicitly the three-loop correlator, which we examined 



in 



]24|1 , in order to understand how we got the formula from the classical solution to 
Heisenberg relation at the (m + l,m) critical point. 

In the second quantized free fermion formalizm ( see, for example, pB| ), 
the connected part of the correlator consisting of the product of arbitrary analytic 
functions f^\A) in two-matrix model can be expressed as 



\i=l 



n . 

n n . 

((N\ II = 4^ : \N)) U / d^K) (X^K) 

1=1 i=l J 

n n 

((N\U--4 i ^--\N))I[W { Hm) , (3-10) 



1=1 1=1 



where 

*(A) = E^a(A), (3.11) 

k 

^(A)= E4^(A) (3.12) 

k 

are second quantized free fermion fields constructed from the orthogonal polynomials 
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£fc and \N)) is a state corresponding to the filled fermi sea, 

a k \N)) =0 for k > N (3.13) 

4 \N)) — for k<N-l. (3.14) 

The normal ordering : ■ ■ • : is with respect to \N)). The connected three-point 
correlator consisting of arbitrary analytic functions f(A), g(A), and h(A) can thus 
be expressed as 

N— 1 oo oo 

(Tr/(i) Trg(A) Trh(A)} = E E E [/(^)U fo(^)]« HA)] U 

i=0 k=Nl=N 
N-l oo AT-1 

E E E[/(^)]*[^)]«[ff(^)]«,(3.i5) 

i=0 fe=iV «=0 

where = (i\f(A)\k). 

Because we are interested in the case of large N limit only, it is convenient to 
use the 'classical' function introduced sect. [2.3.2| . Note that the 'classical' function 
depends on the bare cosmological constant A only through s when we take iV — ► oo 
limit. It is, therefore, legitimate to introduce A(z, s) = lim A(z, s, A) : 

N— >oo 

A(z,8,A) = A(z,s)+0(l/N) . (3.16) 



Since the matrix elements A^ only near the diagonal are not zero, the matrix 
element [f(A)](i)k = [f(A)]i-k,i can be replaced with the coefficient of z k for the 
'classical' function A(z, s = A), 

[f(A)](N) k = ^-J^- 1 f(A(z,s = A)) + 0(l/N) , (3.17) 
at large limit. 

In the right-hand side of eq. ( 13.15 ), the leading terms in 1/JV of the first term 
and those of the second term get cancelled. We have to consider the next leading 
terms. For any integer e, we obtain 

[f(MN + e) k = ±f*Lf(A(z,s = A)) 

df (A(z,A)) 



Ae 1 i dz dA(z,s 



-a 9A 



N 2m J z k+1 ds 
+ (the part independent of e) + 0(1/A^ 2 ) . (3.18) 



The part independent of e comes from the terms 0(1/N) in eq. (|3.17|) . The sec- 
ond term is responsible for the computation in what follows. Using eq. ( |3.18| ) and 
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considering the terms 1/N in eq. ( |3.15| ), we obtain 



Tr/(i) Trg(A) Trh{A) 



d 



^000000 
iV ^ ^ ^ 

<5 1= 0<5 2 =0 5=0 



U52-S 1 ) -[f(A)](N) Sl+ s 2+ i [g(A)](N)ss 2 [/^(AO-^-i 

+(6 + 6 2 + l) [f(A)}(N) Sl+ s 2+1 j- s [g(A)}(N) s _ S2 [h(A)](N) ^ 6l ^ 

+ (5-50 [f(A)](N) Sl+S2+1 [g(A)](N) S -s 2 ^[h(A)](N)^ 5l -i } 
+0(l/N 2 ) . 



(3.19) 



Here we have used the fact that at large N, the original summations in eq. ( |3.15| ) 
can be safely replaced with the triple summations from zero to infinity. 

The summations can be carried out after putting this equation into the form of 
contour integrals using eq. ( |3.17|) . The three-point function eq. ( j3.15| ) in the planar 
limit can thus be expressed in terms of the 'classical' function in the form of contour 
integrals, 

A 1 



Tr f(A) Tig (A) Trh(A) 



N (27T*) 3 i| 21 |>| Z2 |>|23| 



dzidzidz3 



+ 
+ 





(si- 


- z 2 y( Zl - 






Z 2 




(Z2- 


- z x )\z 2 - 


Z3) 2 




Z3 




0*3 




-z 2 y 



6 r 



f(A( Zl )) g{A{z 2 ))h{A{z 3 )) 



f{A{z 1 ))^ K [{A{z 2 ))]h{A{z,)) 
f{A{z 1 ))g{A{z 2 ))^ K [h{A{z,)) 



(3.20) 



where we set A(z) = A(z, s — A). The condition for the contour paths \z\\ > \z 2 \ > 
\z 3 \ follows from the condition that makes the infinite summations converge. 

From the above expression (eq. (|3.20|) ), the three-point resolvent in the planar 
limit is expressed as 



N 
A* 



Tr 



Tr 



1 



Tr 



1 



Pi - A p 2 -A p 3 - A , 



(2iri) 3 J r | 2l |>| 22 |>| 23 | 



d^id^d^ 



Zl 



d 



( Zl - z 2 )\z x - z 3 f OA 
z 2 



1 



Pi - A(zx 
1 d 



p 2 - A(z 2 ) p 3 - A(z 3 ) 



+ 



(z 2 -z i y(z 2 -z 3 y Pl -A(z)dA 

Z3 1 



p 2 - A(z 2 ) 
1 d 



(z 3 - Zl f{z 3 - z 2 f Pl - A( Zl ) P2 - A(z 2 ) dA 



p 3 - A(z 3 ) 
1 

p 3 - A(z 3 ) 



(3.21) 
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where the contour of z\ encircles that of z 2 and similarly the contour of Z2 encircles 
that of z 3 . We pick up only the contributions form the poles at Zi = z* of the parts 
l/[pi — A(zi)]. The contributions from the poles at Zj = z* (i ^ j) give rise to terms 
which do not have the inverse Laplace image. We thus discard these terms. 

We obtain the following rather simple expression for the three-resolvent correla- 
tor, 




94 dz* 3 + ^ + ^ 



a3 d(i d( 2 dCs 1 (zf - z* 2 )(zf - zt) dk x 



(1 <-> 2) + (1 <-> 3) 



A S =A 



(3.22) 



Here, the poles z* are determined through the classical solution to the Heisenberg 
relation and are parametrized as eq. (|3.3|) . 

Eq. ( |3.22| ) agrees with the set of eq. ( |3.2|) and eq. ( |3.7|) , the formula for three- 
resolvent correlator. We have shown how we can get the formula for ra-resolvent 
correlators for the case of n = 3 explicitly. 



3.2.2 Three-loop correlators in terms of loop lengths 



Next, let us consider to how to get the expression for three- loop correlators in 
terms of the loop lengths, performing inverse Laplace transformation with respect 
to the renormalized boundary cosmological constants Q. We will show, later, that 
much physical information can be extracted from the three-loop correlator. The 
generalization to higher loop will be discussed in sect. 4. 



First, we will show we can put eq. ( 3.22|) into a form in which the correlator is 
expressed as a sum of the product of three factors each of which depends only on 
Q corresponding to individual loop. In order to show this, first note that eq. fl3.22j ) 
can be written as 



N / m 1 

T ( Tr 

A *\ px-A 



1 



1 



Tr Tr ■ 

Pi- A P3-A 1 

1 faM\- 2 - l/m Odd 



a 3 2 2 m(m + 1] 



F(9i, 9 2 , 03 



(3.23) 
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where 

1 sinh(m — l)6i 



F(9i, 9 2 , 9 3 ) 



(cosh Q\ — cosh 62) (cosh Q\ — cosh #3) sinh m6\ 
+ (1 <-> 2) + (1 <-> 3) . (3.24) 



Here the following identity is crucial; 

1 / sinh(n — k)a sinh(n — k)(3 \ 

cosh a — cosh (3 \ sinh na sinh n(5 ) 

= _2 yy sinh(w- j-t + l)a sinh(n -j -k + i)(3 
~[ ~[ sinh na sinh n(3 



Making use of the above identity twice, we find that eq. (|3.24j ) is written as a triple 



sum where the summand factorizes into three factors associated with individual 
loops: 

F[8i, 6*2, #3) 

^ m-^ k sinh^ _ tyg 1 s i n h(m — l)9 2 sinh(m — k — j + i)6 3 



k=i j=i i=i 



sinh mQ\ sinh 17162 sinh m9. 



3 



(3.26) 



Here, we should specify the definition of continuum amplitudes at large N. Since 
the leading term of ^Il7=i W(j)i)\ is of order of a~ n n n ~ 2 , where k = a~ 2 ~™ (A*/iV), 
we should renormalize to obtain continuum quantities. The renormalized resolvent 
is defined as 

W ren (Q) = - W( Pl ) = a - Tr , (3.27) 

« « Pi - A 

and the renormalized expectation is defined as 

( ) ren = K 2 ({ )) . (3.28) 

We will omit the superscript ren from now on. The continuum three-loop correlator 
(w + (£i)w + {(-2)w + (£3)) is defined by the inverse Laplace image of the continuum 
resolvent correlator, that is, 



w + {d)w + {C2)w + {Cz) 

dlxj^ d£ 2 J Q ^3e" Cl ^ 1 e- f ^ 2 e-^ 3 (w + (£ 1 )w + (£ 2 )w + (£ 3 )) 
( W + (£> + (£ 2 ) W +(£ 3 ))] , (3.29) 



= £ 
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where w ± {€) is an operator which makes hole with finite boundary (loop) length I. 
Due to the following formula for the inverse Laplace image 



c- 1 



d sinh k9 
d( sinh m9 



sin — KAMI) 

n m 



(3.30) 



where K v (z) is the modified Bessel function, we can obtain the continuum three-loop 
amplitude in a rather compact form: 

w + (e 1 )w + (£ 2 )w + (£ 3 ) 



1 



-2— ±- 



m(m + 1) V 2 / 

(A/f\3 m—lm—k k 

— ) hhh Y] Y) £ Km-k(M^) K m-j-t+i (M£ 2 ) K m-k-i+i (Mi*) 
2 ' fe=l j=l i=l 



where we introduced a notation, 



K r( Mt) . K r (Ml) 

71/2 



(3.31) 



(3.32) 



The expression for the summation in eq. ( |3.31| ) looks asymmetric with respect to 
the loop indices. By elementary algebras, we can convert it into a form which have 
explicit symmetry with respect to the interchange of loops: 

w + {h)w + {t 2 )w + {h) 



m(m + 1) V 2 / 



MV 



-1*2*3 



£ K 1 _k±(M£ 1 ) K x _ k _,{Mt 2 ) K^Mi, 



( fel _l,fc 2 _l, fc3 _l) 



(3.33) 



Here we have denoted by T>^ 



V 



(m) 



< (a 1; a 2 , a 3 ) ^ — a j > for j = 1 ~ 3 

I i&i) 
3 

2_, Q-i = even < 2(m — 2) , aj = 0, 1, 2, • • • 



(3.34) 



Eqs. (|3.33| ) and (|3.34|) give the final expression for the three loop correlator at 
the (m + l,m) critical point. It is interesting that the selection rule obtained in 
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eqs. Q3.33] ) and ( |3.34j ) agrees exactly with the fusion rules for the diagonal primary 



fields in the Kac table of underlying conformal field theory of the unitary minimal 
model (m + l,m) |L7[]. In fact, the fusion rules for the diagonal primary fields read 
as 

(4>a 4>jj 4>kk) ^ , (3.35) 
if and only if i + j > k + 1 and two other permutations and i + j + k (= odd) 

r 



< 2m — 1 hold. This set of rules is nothing but V^ mj 



3.2.3 Boundary conditions of loops 



A similar expression to eq. ( |3.33| ) was obtained in ||34|| , where loop correlators were 



examined for closed string with one- dimensional discrete target space, that is, the 
degrees of freedom for matter part are labeled by the points of Dynkin diagram. The 
matter degrees of freedom are labeled also by the discrete momentum p instead of 
the discrete target space coordinate x. They examined the loop correlators treating 
the boundary condition of each loop £j to be specified by a single momentum pi. 
Thus, it follows directly that the three-loop correlator which is specified by three 
momenta pi,p 2 and p% and loop lengths, is proportional to the expectation value for 
wave functions for matter part, that is, 

CpiP2P3 = ^ S(pi)S(p 2 )S( P3 )/ S x > (3.36) 

X 

where 5? % and S x are the wave function of a point particle moving on the discrete 
target space with momentum p^ and that for ground state respectively, 

/ 2 \ 1/2 

= ( x \p) = (——A sinvrpx (3.37) 

x — 1, 2, • • • , h — 1 

p= l/h,2/h,---,{h-l)/h . 

For example, for the Ah-i Dynkin diagram, each momentum takes discrete values 
of r, ■ ■ • , and C PlP2P3 is nonvanishing only when hpi satisfy the equivalent rule 
we found in eq. ( |3.33| ) and eq. ( |3.34j ), that is, 



(h Pl - 1, hp 2 - 1, hp 3 — 1) G Vi h) . (3.38) 

The similarity between the three-loop correlator in the case of the closed strings 
in discrete target space and that in the two-matrix model we found in eqs. Q3.33j ) 
and ( |3.34j ) indicates that the each terms in the sum in eq. ( |3.33| ) represents the 



24 



amplitude with the the loops specified by the momentum ^. It appears that we can 
decompose the loop operator in the two-matrix into parts each of which specified 
by a momentum — : 

J m. 



m—l 



u w + (£) ~ £ 4Mty , (3-39) 



k=l 



From the selection rules in eq. fl3.34|) , we can deduce the fusion rules for the 



gravitational descendants as well as for the gravitational primaries. In other words, 
some fusion rules are satisfied among all of the scaling operators including the grav- 
itational descendants as well. We suggested first that the selection rules in the 
three-loop correlator correspond to those for the gravitational primaries by exam- 
ining the limit of small loop length |23fl. In ||26|| , we obtained the fusion rules for 
all of the scaling operators from the three-loop correlator. We will discuss these 
issues in detail in later subsections after we examine the two-loop amplitude and an 
expansion of the loop operator in terms of local operators. 

3.3 Expansion of loop operators 



In [[30], |31| it was proposed that a loop operator can be replaced by a sum of local 
operators if the loop correlator does not diverge when the loop shrink to a point. 
This was discussed explicitly in the case of the one- matrix model and c = 1 case. 
We apply this idea to the general minimal models coupled to gravity and discuss 
the correlation functions for the scaling operators. In order to derive the form of 
the expansion of the loop in terms of the local operators, let us first consider the 
two- and one-loop amplitudes. 

3.3.1 Two-loop correlators from three-loop correlators 

Let us derive two- and one-loop correlators from the three loop correlator eq. (|3.33| ) . 
As we shrink one of the three loops, the three-loop correlator should approach the 
derivative of the two-loop correlator with respect to the cosmological constant. Con- 
sider shrinking the third loop M£ 3 in eq. ( |3.33|) . Since, for Mi < 1, we have 



/ m / v m m J 




I_ l+]L {M£) 

(3.40) 
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the leading contribution in the summation in eq. ( |3.33| ) comes from the part of 
&3 = 1 and we have 

i fm m—1 

« 7 -7SWT/ K x _±{Mi x ) K x _±{Mh) , (3.41) 

m[m + 1) L (1 / m) r^. ™ «> 

for M£ 3 <C 1. This should be proportional to the derivative of two- loop correlator 
with respect to \i. In fact, by the explicit calculation similar to the case of three-loop, 



one can obtain |15|, |23 

d 



-y m—1 



m[m + 



rr^E K^Mh) K x _u{Ml 2 ) . (3.42) 

1J fe=1 



It is clear that eq. ( |3.41| ) and eq. ( |3.42| ) are consistent. 



Due to a relation ^- = M q M and a formula of a integral (for a ^ 

J dzzK u (az)K u ((3z) = ^ r ^—^{f3K v (az)K u _ l (f3z) - aK v _ l (az)K v (fiz)}, 

(3.43) 

one can obtain the two-loop correlator (for l\ ^ £ 2 ): 



w + {£ 1 )w ± {£ 2 ) 



1 M 



m—1 



m 2 l\ - t 



— / ~ V m m 



2 fe=l 



1 M /,/, 
mY£ 1 + (±) m £ 2 



in m J 

E (±) fc_1 #±(M4) K 1 _k(M£ 2 ) . (3.44) 



fc=i 



3.3.2 One-loop amplitudes from two-loop correlators 

Shrinking l 2 in eq. (]3.44|) as well, one should have the derivative of the one-loop 



amplitude with respect to the cosmological constant. For Mi 2 < 1, we have 

JL i 

. 1 Dm / fjf\ „ . 

w +(*>±(* 2 )) « - K±(M£i) . (3.45) 
' ml(l/m) \ 2 / ™ 

In fact by explicit calculation, one can obtain 

= ™ (t) ^ (m «- (3 - 46) 
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Note that this amplitude is nothing but the wave function of the dressed identity 
operator. Performing the integral with respect to /i, one obtain the one-loop ampli- 
tude: 



1 

1 + - 

m 



1 +— 

_ x sin7r/m ^ 



ir/2 V 2 



1+: 



.{Ml x ) 



(3.47) 



Note that l\ (w ± (£i)) is the wave function of a boundary operator di +m |33|] which 
couples to the boundary of two-dimensional surfaces (i.e. loops), 



(M 



1 + ™At> K ^ m - 



(3.48) 



3.3.3 Expansion of loops in local operators 

In [HO], in the case of one- matrix, it was argued that the loop operator can be ex- 
panded in terms of local operators inside the loop correlators, that is, the loop can be 
replaced with the infinite combination of local operators, except some special cases. 
Whether this replacement can be done safely or not is connected with whether the 
corresponding classical solution exits or not in the limit of small length of corre- 
sponding loop. This claim is quite natural because, in the one-matrix model, all of 
the scaling operators are expressed in term of one matrix $ as 

(Tj = Tr(l - <t>y +1 / 2 = ]T a n (j)n- 1 Tr $ n . (3.49) 

n 

On the other hand, in the two-matrix model, since there are two kinds of matrix A 
and B, there can exist many kinds of microscopic loops, Tr(A ni B n2 A n:i ■ ■ •). In the 
case of the two-matrix model, thus, the direct connection of the scaling operators 
to the loops Tr A l l a or Tr B e '/ a is not clear. At first sight, the expansion of loops 
in terms of local operators is not legitimate. We think, however, this expansion 
is possible by the following reason. When one of the loops on two-dimensional 
surface shrunk to a microscopic loop, the loop represents local deformation of the 
surface. The microscopic loop can be considered to be replaced by the insertions of 
local operators. As we will see later, the loop correlators except one-loop case are 
continuous when the length of one of the loops approaches zero, so that we expect 
that a macroscopic loop can also be replaced by a sum of local operators. 

In the following, we find the form of the expansion of the loop operator. First, let 
us represent the two-loop correlator in term of off shell renormalizable wave function 
§0, 

(£sinhvr£) 1/2 K iE {M£) . (3.50) 
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We show that the two-loop correlators are expressed as 



00 , „ E sinh it E 



/ dE— — F K iE (M£ 1 ) K iE (M£ 2 ) (3.51) 

Jo cosh T\E — COS 7T — 



From eq. ( |3.42j ), it is reasonable to assume 



m—1 „• _ k 



= J] (i)^ 1 ^!^ Wfc (£)" (3.52) 



sin 71" — 



k=l 



tt/2 



in the case of two- and three-loop correlators, where we have introduced loop op- 
erators Wk(£) which represent loops with some distinct matter boundary condition 
(see sect. |3.2.3|) . From eq. ( |3.42|) , we have 



"' 1 / sin ;•,-- * ' 



(w+MwH**)) = E (±)^ ( ~^jf ) (MUM**)) (3-53) 



k=l 

and 



° {w k {h)w k {i 2 )) = --^-£i£ 2 K x _±{Mt x ) K x _u.{M£ 2 ) . (3.54) 



dM x v ' v " m 2 
Making use of a formula 

and replacing t with tM 2 , we have 

M£ x l 2 K x _k_{Ml\) K 1 _±(M£ 2 ) 

m m 

Carrying out the integral with respect to M, and from eq. Q3.54|) , we have 

H(^iH(4)> = -L 
4m 



Due to a formula, 



/■°° , £sinli7ri? . , . 

zK l _\ p \{z)= / dE K iE {z) , 3.58 

jo cosh ttE — cos 7rp 
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the right hand side of eq. ( |3.57| ) turns into 

1 r°° dt 
4m Jo t 

E sinh ttE 



dE- 

i ) cosh -kE — cos np 



Ki 



-1*2 



tM 2 



exp 



i / 1 \ 2 

Using a formula eq. (|3.55| ) again, eq. (|3.59|) turns out to be 



2t 



1 r 00 
— / d£- 



sinh 7ri? 



} ._ ,. — K iE (M£i) K iE (M£ 2 ) 

2m Jo cosh 7r_E — cos 7r— 



(3.59) 



(3.60) 



Putting eq. ( |3.60|) and eq. (|3.53|) together, we have proved eq. (|3.51 



Let us go back to eq. ( |3 . 5 1| ) and perform the i?-integral. The integral can be 
carried out by deforming the contour. The residues for poles 



E = ±i{— + 2n) 



n 



0,±1,±2,--- 



m 



(3.61) 



contribute to the integral and, after all, we obtain the following expansion for the 
two-loop correlators (for l\ < £ 2 ) 



{w k {i x )w k {l 2 )) 

i oo 

1 E ( 



k \ /sni7r— 

— + 2n ) 

m > \ it 2 



(3.62) 



and 



m—l oo 



m 



k=l n=— oo 



k 

— + 2n 

m 



Y2n 



(Mi,) K ±+2n (M£ 2 ) . (3.63) 



Since the two-loop correlators eq. ( |3.44|) or eq. ( |3.63| ) do not diverge when the 
length of one of the loops approaches zero, we expect that one of the loops can be 
replaced by an infinite combination of local operators. From the consideration of the 
minisuperspace Wheeler-deWitt equation and the scaling behavior, we expect that 
the wave function of the scaling operator dj is proportional to M^ m Kj/ m (M£). 

indicate that the following normalization of the wave 



Eq. and eq. 

function 

w + (£) 



J_ 

m 

J_ 

m 



M\ ™ sin — 7r , 

1 / 7T/Z m 

— K±{M£) , j > 1 ^0 (mod m) 

2 J m 



(3.64) 
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would be reasonable. Note that the normalization factor sin ^ir in eq. ( |3.64| ) is 
consistent because there are no scaling operators dj for j = (mod m) in the matrix 
model. Comparing eq. (|3.64j) with eq. (|3.63|) , we expect the following expansions of 
the loop operators in term of the local operators: 

-r m— 1 oo , tl/t v — \—+2n\ 

" W± W = ~ E E (±)* _1 (-5-) " ^ +2n ,(M£) a, fe+2mn | " . (3.65) 

171 k=l n=-oo W / 

These expansions are the generalizations of those in the case of the one-matrix model 



30 1 to the case of arbitrary unitary minimal model coupled to two-dimensional 
gravity. 

Since the loop correlators are symmetric under the interchange of two kinds of 
loops, that is, (w + (£i) w + (£ 2 )) = ( w ~(h) w ~(h)), the wave functions of the scaling 
operators with respect to loop w~(£) are read as 

(a 3 w-(£)) = (-iy- 1 (a 3 w + (£)) . (3.66) 

The wave functions with respect to the loop Wk{£) are 

(k \ /M\ l - +2n| 

(9\ k +2mn\Wk'(i)) = <W [- + 2n j %+2n|( Mf ) " ( 3 ' 67 ) 



3.4 Relation to the mult i- matrix model 



Let us comment on the relation between loops in the two-matrix model and those 
in the multi- matrix chain models. The lowest critical point of the (m — l)-matrix 
chain model represents also the (m + l,m) minimal model, which corresponds to 
A m _i Dynkin diagram, coupled to two-dimensional gravity. From the observation, 

m— 1 ,,~ 

l* = l> = ElP^lfc) amtfi . (3-68) 

h=l 
m—1 

\x = m-l) = Eb=^)(^i) sinTr^(m-l) 
fc=i 

m-l 

= Eb4)fe) ("l^sinvrA , (3.69) 

fe=i 

where (x\p) is the wave function introduced in eq. ( |3.38|) , we think w + (£) and w~(£) 
should correspond \x = 1) and \x — m — 1) respectively. The loop operator w^(£) 
created by the x-th. matrix of the (m — 1) -matrix chain model, thus, corresponds 
to 

m— 1 -.,„ 

W=EIp4)(^) , (3.70) 

k=i 
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and would be represented accordingly as 

m ~i sin7T-— "X 

w {x) (£)=T, TK-Mf) , x = l,---,m-l . (3.71) 

k=i n / 2 

We think this relation is valid at least for loop correlators with less than four loops. 
Using the relation, we can construct the loop correlators of the multi-matrix models 
from those of the two-matrix model. 



3.5 Three-point functions and fusion rules 

3.5.1 One- and two-point functions 

Let us consider the correlators of the scaling operators. We can extract these from 
loop correlators due to the relation eq. (|3.65| ). 



Since the one-loop amplitude diverges when the loop length approaches to zero, 
this amplitude include the contribution which is not represented by the local oper- 
ators. Putting the one-loop amplitude into 



= fv) 2+ " 1 ( J 2+i( M ^) " !-2-±(M£) - n(M£) + I_ 

\ £ J \ m m m 



(3.72) 



and extracting the parts proportional to l v {y > 0), which parts would be considered 
as the contributions from local operators, we can obtain the one-point functions of 
the scaling operators: 



= > ( 3 - 73 ) 

(a 1+2m ) = mi^—j , (3.74) 
= 0, j^l,l + 2m. (3.75) 



Let us turn to the two-point functions. Substituting eq. ( |3.65| ) into eq. ( |3.64j ), 
we obtain the two-point functions: 



/M\ 2j/m 

I — \ , «,j^0 (mod m) 



{didj) = 6ij j 
Note that we obtain diagonal two-point functions. 



(3.76) 
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3.5.2 Three-point functions 



As for three-point functions, using the formula, 

E sinh ttE 



zKx-\ p \(z) = / 

•J 



dE r K iE {z) 

o cosh tiE — cos 7rp 

oo 



- E 

n=— oo 



|p + 2n| 



I i o i-|P+2n|l 

sm7r|p + in I 

we first expand the three-loop correlator eq. ( |3.33| ) as 

7 m(m + 1) V 2 J 
k ± + 2ni) (J + 2n 2 ) (| + 2n 3 ) J,^ +2ri (M*)!^ (M^)!*. 



(3.77) 



E E E E 

X>3 ni=-oo ti2=-oo «3-oo 



■2n 3 ] 



(M4) • 
(3.78) 



Comparing eq. (|3.78|) with eq. (|3.65 ), we can extract the three-point functions 



\ cr |fci+2mni|C"|fc2 +2mn 2 1 (J \ k 3 +2mn 3 1 



m(m + 1) 



n(^ + 2mn,)( — ) , (3.79) 



8=1 



2 / 



where 



Cfclfc2&3 



1 





(fel - 1, fc 2 - 1, &3 



1) G Pf } 



(3.80) 



otherwise 

For rii — 0, eq. (|3.79 ) is nothing but the correlator of the gravitational primaries. 
For the gravitational primaries, eq. (|3.76|) and eq. Q3.79|) agree with the correlators 



obtained in [£2| from the generalized KdV flow up to a factor —2. Note that we 
obtain, here, the correlators of the gravitational descendants as well. 



In ||19|| , the fusion rules for the gravitational primaries were examined in contin- 
uum framework. Note that we have found here the fusion rules for the gravitational 
descendants as well as for the gravitational primaries. These fusion rules are sim- 
ilar to those for the gravitational primaries due to the factor C^^ka in eq. ( |3.79| ). 
Introducing the equivalence classes by the equivalence relation 

Sfc ~ ?|fe+2mn| , n e 2Z , (3.81) 

we can consider the fusion rules in eq. fl3.79|) as fusion rules among [a^] (k = 
l,---,m — 1). Note, here, that the class does not correspond to the set 
which consist of the gravitational primary and its gravitational descendants 
ai{O k ) ,1 = 1,2, 



m 



p2fl introduced from the viewpoint of KdV flow. The three- 
loop correlator eq. (|3.33j ) represents the fusion rules for all of the scaling operators 
including the gravitational descendants in a compact form. 
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3.6 Further on the fusion rules 



Let us examine the fusion rules in eq. ( |3.79| ) further and consider the relation of the 
scaling operators to the primary fields in the corresponding conformal field theory. 

In the (p, q) minimal conformal model, the primary field $ rs has the conformal 
dimension 

A r ,= ipr ~ qS? - (P - 9)2 , (3.82) 
Apq 

where r and s are positive integers. Since 

^r,s ^r+q,s+p ^q—r,p—s j (3.83) 

the corresponding primary fields can be regarded as the same one. The integers r 
and s can thus be restricted in the range 

1 < r < q — 1 

1 < s ^ (mod p) (3.84) 
pr < qs 

(see fig. [| ). In fig. [^, the primary fields in the region ((2)) or ((2))' are the secondary 
fields of those in the region ((1)). In general, the fields in the region ((n + 1)) or 
((n + 1))' are the secondaries of the fields in ((n)) or ((n))'. Since the secondary 
fields correspond to null vectors, those fields decouple. One can thus construct 
consistent conformal field theory which include only the primary fields in the region 
((1)) (i.e. inside the minimal table), that is, the minimal model (p,q) fL7| . Coupled 
to Liouville theory, however, the fields outside the the minimal table fail to decouple 
]I8"| and infinite physical states emerge accordingly. These states are considered 
to correspond to the primaries outside the minimal table. This correspondence is 
implied by the BRST cohomology [BO], of the coupled system. 

Denoting the gravitational dimension of the dressed operator for $ r s by A^ s = 
1 — in the minimal model coupled to Liouville theory, the following relation was 



shown 20, 21 



a r ,s p + q- \pr - qs\ 



(3.85) 
7 2g 

where r and s take the values in the range eq. (|3.84j ). On the other hand, in the 



matrix model, the corresponding relation for the scaling operator dj is 

otj _p + q-j 



(3.86) 



7 2g 

From eq. ( |3.85| ) and eq. Q3.86j ), we should take as 

j = \pr — qs\ , j = 1, 2, • • • ^ (mod q) , (3.87) 
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for (7j. 

Consider now the relation of d\k+2mn\ to the primary field $ rjS of the unitary 
(m + 1, m) minimal model. Let us first compare the two sets 

S k = [\k + 2nm\ \ neTl] , (3.88) 

and 

| \pr — qs\ | = | |(m + l)r — ms\ j 

= { r'+ (a -r - l)m } , (3.89) 

where r and s are positive integers in the range 



(3.90) 



and r' = m — r. Note that we include s = (mod m + 1) here. Dissolving the set 
Sk into two sets as 

S k = St<BS£ , (3.91) 

where 

S£ = {k + 2nm \ n = 0,1,2,- ••} , 

= {(m -fc) + (2n' + l)m | n' = 0,1,2,- ••} , (3.92) 

and comparing eq. ( |3.92|) and eq. ( p.89|) , we can express the sets S£ and in terms 
of \{m + l)r — ms\ as 

S k = | |(m + l)r — ms\ \ r' — k, s — r — 2n + 1, n = 0, 1, • • • j , 

S k = { \(m + l)r-ms\ \ r = k, s-r = 2n' + 2, n' = 0,1,- ••} . (3.93) 




From eq. (|3.93| ), the following correspondence is obtained: 

&\k+2mn\ (n>0) $ m -fe,r+2n+l > 0) 

V\k+2m(-l-n')\ ( n> > 0) ^ $fe, r+2n'+2 («' > 0) , (3.94) 

where s 7^ (mod m + 1). 

In [j33|l , it was suggested that the scaling operators a^, j = (mod m + 1), 
should be identified as the boundary operators which couple to the boundaries of 
two-dimensional surface. These scaling operators do not have their counterparts in 
the BRST cohomology of the system coupled to Liouville theory. From 

\(m + l)r — ms\ = (s — r)(m + 1) — s = (mod m + 1) , (3.95) 
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r 




Figure 2: the range of (r, s) 



s = (mod m + 1) in the range eq. (|3.90|) should correspond to the boundary 
operators. 

As an example, we depicted the scaling operators on the r-s plane for the case 
of m = 4 in fig. [3[ 

Concerning the scaling operators inside the minimal table, the fusion rules in- 
volving non-diagonal operators dose not agree with the fusion rules of the unitary 
minimal model; 

($ ri)Sl $r 2)S2 ^r3, S3 ) ^ iff 

(n - l,r 2 - l,r 3 — 1) 6 V { t ] and {s x - 1, s 2 - l,s 3 - 1) G (3.96) 

For example, in the Ising model (m = 3) the three point function for the energy 
operators vanish, 

($i,3$i,3$ 1,3) = • (3.97) 
Coupled to gravity, however, the corresponding three point function does not vanish. 



3.7 Boundary operators 

3.7.1 Boundary operators and touching of loops 



In it was proposed that the scaling operators which do not occur in the BRST 
cohomology of Liouville theory are boundary operators and one of them, which is 
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Figure 3: the scaling operators a\ k+2m n\ m the (5,4) minimal model coupled to 2D 
gravity 

S3 = diiOi) in the case of pure gravity, was in fact proven to be a boundary operator 
for the one-matrix model and the Ising model case. We would like to examine the 
role of the operators a n ( m +i), n = 1, 2, • • • 7^ (mod m) as well as cx m +i for general 
unitary minimal models. 

Let us denote these operators by 

Bn = Vn( m +i), n = 1, 2, • • • ^ (mod m) . (3.98) 

In the matrix models the loop amplitudes contain the contribution from the con- 
figuration with loops touching each other. In two-loop case, let us consider the 
configuration in which the two loops touch each other on n points. When we shrink 
one of the loops to a microscopic loop, the other loop splits into n loops, which 
are stuck each other through the microscopic loop (see figs.|], [| and ^). Since the 
microscopic loop represents a sum of the scaling operators, the wave functions of 
some scaling operators contain the contribution from the surfaces with split loop. 

We now show that the boundary operators indeed represent these configurations. 
From eqs. (|3.64|) and ( |3.47|) , the wave function of B n and the one-loop amplitude 



are 

[B n w + (£)) = n(l + -) ( T J K n(i + ±)(Me) , (3.99) 

(w+w) = (i + l) r 1 (f ) 1+ 'k 1+ ±(mi) . (s.ioo) 



36 



In the space of Laplace transformed coordinates, we have 



£ 



2 J 



2 cosh n(m + 1)9 



£ 



M 



— ) 2cosh(m+l)6 l 



where we have used the relation 



£ 



-t- x \v\K v {Ml) 



2 cosh mvQ . 



(3.101) 
(3.102) 

(3.103) 



Note here that w + {€) represents a loop with a marked point and £~ 1 w + (£) repre- 
sents a loop without a marked point. Since cosh n(m + 1)0 can be expressed as a 
polynomial of cosh(m + 1)6, 



2 cosh n(m + 1)0 



,(») 



2 T n (cosh(m + 1)0 

[n/2] 

^ C ( n) [2cosh(m + 1)0 

r=0 

(— l) r n /n — r 
n — r V r 



n— 2r 



(3.104) 



where T n is the Chebeyshev polynomial, we obtain the following relation: 



E 

r=0 



2 / 



{£ [-(w+ (*)>]}' 



n-2r 



(3.105) 



In the space of loop lengths, the above relation means that the wave function of B n 
is equivalent to a sum of the convolutions of disk amplitudes: 



B n w + {i) 



E 

r=0 



2 / 



-1 



in— 2r 



n— 2r 



(£) . (3.106) 



Here we introduced a notation ^4^ = (u» + (£)), and OAf (£) denotes the convolu- 
tion of s Af(£Ys, for example 



o Jo 



M X M 2 5{h + £ 2 - £)At(£i)At(i2) 



(3.107) 



From eq. ( |3.106| ) we can conclude that the operator B n couple to the point to which 
s (< n) parts of the loop are stuck each other in the case of one-loop amplitude. 
When there are more than one loop, we infer that the operator couples to the point 
to which s parts of several loops are stuck each other; the operator will not recognize 
that it is touching different loops this time. 
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Using the following relation 

[(n-l)/2] 



E 



r=0 



2 cosh x 
we also obtain 

i\QAt] n (£) = (-1) 



?? 



2 cosh(n — 2r)x, (up to constant) 



[(n-l)/2] , . M 2r(l+i) 

_ ir+ i E 



r=0 



r/ V 2 J 



B n „ 2r w + (£)) . 



(3.108) 



(3.109) 



Here we drop the constant term in eq. (|3.108|) when we carry out the inverse Laplace 
transformation. From eq. (|3.109| ), we see that the boundary operator coupled to the 
point on which n parts of loops are touching each other is given by 



[(n-l)/2] 
. 1)B+ 1 £ 

r=0 



2 J 



B 



n—2r 



(3.110) 



Now let us consider the boundary operators when there are two loops on two- 
dimensional surface. As for Bi, we expect that (w + (£\)w + (£ 2 )£>i) should be propor- 
tional to (£i + £ 2 ) {w + (£\)w + (£ 2 )) . Let us confirm this in the following. 

From the three loop correlator ( |3.33| ), the expansion of loop operator ( |3.65|) and 
the wave function of &\k+2mn\ (|3.64j) , we obtain the following correlator with two 
loops and a local operator: 



W + (£ 1 )w + (£ 2 p\k 3 +2mn 3 \ 



m + 



k\,k2 



x£ x £ 2 {— + 2n 3 )K lhL {M£ l )K 1 _ tz {M£ 2 ) 

7YI m m 



(3.111) 



Consider the amplitude for B\ = B\ = <r m +i = <Xim=i_ 2 |. Since Ck l k 2 m-i is nonvan- 
ishing only for the case of k\ + k 2 = m, we obtain the following amplitude: 



m— 1 



w + {£ 1 )w + {£ 2 )B 1 ) = - hr ^K^M^K^^M^) . 

\ 2 J m m 



m 



(3.112) 



Comparing eqs. ( |3.112| ) to ( |3.44j ) we obtain the desired relation: 
V(4)w + (4)£i) = (4 + 4) (w + {£x)w + {£ 2 ) 



(3.113) 

Next, let us consider B 2 . Since we infer that the insertion of B 2 should play the 
role of connecting two parts of loops together, we expect the following relation: 

w + (4)^ + (4)£ 2 ) = 24 

«2 



+2£ 2 



d£[ {w + (£' 1 )w + (£ 2 )){w+{£ 1 -£' 1 ) 
d£! 2 (tw + (4)tw + (0) ( w+ (4 - Q 



+2£ 1 £ 2 (w+(£ 1 +£ 2 )) . 



(3.114) 
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Figure 4: (a): A surface with two loops touching each other on a point, (b): When 
one of the loops shrinks to a microscopic loop the microscopic loop is equivalent to 
the insertion of the operator denoted by the dot on the loop. 




(a) (b) 

Figure 6: The case of a surface with two loops touching each other on three points. 
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The third term in the right hand side of eq. ( |3.114j ) represents the contribution from 



the surfaces with loops w + (£i) and w + (£ 2 ) touching with each other on a point. 
Let us confirm the relation (|3.114j) in the following. In this case, it is convenient to 



consider in the space of Laplace transformed coordinates Q. In this space eq. ( p.lll| ) 
reads as 

^+(G)^+(C 2 ) a |fc3+2mn3 |) = (y) Q + 2 ^3) 

_d d_ J ^ sinh(m -jijgi sinh(m - k 2 )9 2 \ , , 

X ^Ci5C 2 \^ 2 lk2k3 sinhm^ sinhm# 2 J" [6 ^ 

Due to the relation 

sinh(m — kx)9i sinh(m — k 2 )9 2 



, , sinhm^i sinhm#2 



— 1 / sinh(m — k%)9i sinh(m — k^)9 2 



2(cosh^i — cosh #2) \ sinhm^i sinhm#2 

we have 



(3.116) 



1 /AA^- 2+ l^+ 2 " 3 l ,k 



1 2{m -hi) V 2 / m 
d d j 1 / sinh(m — k%)6i sinh(m — k 3 )9 2 \ 1 

dC,i ^C,2 \ cosh 61 — cosh 9 2 \ sinh mB\ sinh m6 2 ) J 

(3.117) 

Since we should take k 3 = 2 for B 2 = — &2(m+i) (for m > 3), we obtain the amplitude 
for B 2 



^ + (Ci)^ + (c 2 )s 2 ; 

-1 /M\» d d \ 1 /sinh(m-2)0i sinh(m - 2)# 2 



m V 2 / d(i d( 2 [ cosh 9\ — cosh 9 2 \ sinh m9\ sinh m# 2 y J 

(3.118) 

On the other hand, from the amplitudes 



^ + (Ci)^ + (C 2 )) = ^m -sh^-cosh^ 



<9Ci 9(2 cosh m9\ — cosh m9 2 
d \ 1 sinh^x 1 



dC,2 1 cosh 9\ — cosh 9 2 mM sinh m9\ Ci — C 2 



(3.119) 



iy+(C)) = -(— 2cosh(m + l)#, (3.120) 



2 y 
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we obtain the following relation 
-2- d 9 



{(^ + (Ci)^ + (C 2 ))(^ + (Ci))} 

^ + (d)) - (^ + (c 2 ): 



2) 



C1-C2 



2 /M\~ 9 d 



I sinh 6 1 ! cosh(m + 1)9 \ 



m \ 2 / 9Ci <9£ 2 [ cosh 9\ — cosh #2 \ 



sinh m9\ 



(1 




One can easily show that the right hand side of eq. ( |3.121| ) agrees with that of 
eq. fl3.118p . Putting eqs. (|3.121|) and (|3.118| ) together and performing the inverse 
Laplace transformation, we obtain the desired relation eq. ( |3.114j ). 

We have shown that the operator B 2 connects two parts of loops together in the 
case with two loops. We infer that similar phenomena occur in general; the operator 
B n would connect n parts of loops together in the case with any number of loops. 



3.7.2 Connection to the Schwinger-Dyson equations 

We can observe close relationship between the boundary operators and the Schwinger- 
Dyson equations proposed in [|K|. Continuum limit of the Schwinger-Dyson equa- 
tions for loops in the two- and multi-matrix models were proposed in [O] under 



some assumptions. It was shown [K| that these equations contain W3 constraints, 
which were derived explicitly in . The integrability of these equations were shown 
in [|nj. These facts justify the proposed Schwinger-Dyson equations. 

Let us consider the connection of the boundary operators with the Schwinger- 
Dyson equations. For the (m+1, m) minimal models, the following Schwinger-Dyson 
equations were proposed in fllQ| : 

j*dl' (w^(i')w^ (£ - £'; [H(a)] j ) • • 

+gJ2ii (i + U\ [H(*)Y) w^(h) ■ ■ ■ u)W(4.i)tfl«(4fi) " " " w (1 \Q)' 

i 

+ (w® (t, [H(a)] j+1 ) w {1 \l x ) ■ ■■w {1 \l n ))' « , 

for j = 0,---,m-2 , (3.122) 

and 

(w {1} {£; [Hia)}™- 1 ) w^\i x ) ■ -/'(y)' « . (3.123) 
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Here (• • •) represent loop correlators that are not necessarily connected, and 
w^ l >{£) represents a loop operator corresponding to a loop created by the matrix 
in the multi- matrix model. The operator 7i(a) describes an operator which 
changes the 'spin' on a loop locally from 1 to 2. Also tyW (£; ['H(cr)p') describes a 
loop with [H(a)Y inserted. The symbol ~ means that as a function of £, the quantity 
has its support at £ = 0. 

From eq. ( |3.122| ) for j = and n = 1, we have the relation 

h (w w (4; H(a)) » (1) (l 2 ))' + £ 2 lw^\h)w^ (£ 2 ; H{a)))' 

+z l p d£\ (w^{e[) w ^{e 1 - e x )w^\i 2 ))' 

+£ 2 fj d£' 2 (w^{£ 1 )w^(£' 2 )w^{i 2 - £' 2 )) 

+2g£ 1 £ 2 (w^(£ 1 + £ 2 ))' ^0 . (3.124) 

The planar part of the above relation agrees with eq. (|3.114 ). Note that the loop 
amplitudes in eq. ( 3.114j) represent connected correlators. 

This agreement implies that TC would correspond to B 2 . Taking into account the 
fact that B n (n = mod m) do not exist and eq. ( |3.123| ), it is legitimate to consider 
that the amplitude (for j = 1, • • • , m) 

{w + {£ l )---w + {£ n )B ] ) (3.125) 

corresponds to the connected part of the amplitude 



J2 fda t (u^ft) ■ ■■w^(£ i ^)w^ (4 Mai)]'- 1 ) w^(£ i+1 ) ■ ■■w^{£ n ))' . 

(3.126) 



i=l 



42 



4 Multi-loop correlators 



In this section, we generalize the discussion in sect. ^ to the cases of higher- 
loop. First, we derive the formula of the n-resolvent correlators, which we quoted in 
sect. pTT| , and point out that the structure corresponding to the crossing symmetry 
of the underlying conformal field theory can be seen in the loop correlators. We 
then discuss the four-loop correlator in detail. 

4.1 The derivation of the n-resolvent correlators 

Consider in the two-matrix model the connected part of the correlator consisting of 
the product of n-resolvents. 

Ir^I ■ <4i) 

/ -. \ n— 2 

It should be noted that this expression is at most ( J due to the large N factor- 
ization of the correlator consisting of the product of singlet operators. In the second 
quantized notation, eq. (|4.1|) is expressible as 



n r 1 
((N\U: d^\\ i )-—*(\ i ):\N)) 

i=l Pi ^i 

n n f ~ \ 

= ((N\ II ■ 4,% : |JV» II / d^ ki (\i)-— T t ki (\ i ) 

8=1 i=l J Pi- A « 

n n 1 

= «n\ n = : m n^ir— > (4 - 2) 

i=l i=l Pi A i 

The normal ordering : ... : is with respect to the filled sea \N)). We introduce a 
notation 

— ^- ( Zi ;A i ,A,N) = J2 z i<3i-t\-^\3i> (4-3) 
p — A\ g p — A 

A l =j i A/N = A + A~ Jl /N . (4.4) 

The evaluation of ({N\ n"=i : <A, a k ■ 1-^0) by the Wick theorem provides (n — 1)! 
terms of the following structure: each term is given by the product of n-Kronecker 
delta's multiplied both by a sign factor and by the product of n-step functions to en- 
sure that the summations over the n-indices j\,j2 ■ • • and j n are bounded either from 
below (> 0) or from above(< —1). We denote this product by 0(ji,j2,- ■•jn' 1 &). 
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These (n — 1)! terms are in one-to-one correspondence with the circular permuta- 
tions of n integers 1, • • • , n, which we denote by S n . The a is an element of S n . For 
large N, we find 

E^B^,-i,», w (n/|)ni(?)'" 



X 



1 / n ~ 8 \ n 2 n 1 

n K ,-X(„,;A„) Iv-A-^dW • (4-6) 

Note that in the large N limit, we can use — — ^ in place of — ^ (zi, A i? A, N) 
according to the same reason as stated in the case of three-loop correlator. The 
sgn(a) denotes the signature associated with the permutation a. 

Let us define 

m\D m (z,z') = -J2r(z'/z) ] £ r{z'/z) j , m = 0,.-. . (4.7) 

In the continuum limit we will be focusing from now on, it is sufficient to use 

D m(z,z') « {z _ l zl)m+1 = D m (z - z>) . (4.8) 

Let sgrii(a) be +1 or —1, depending upon whether the restriction on the summation 
over ji is bounded from below or from above respectively. It is not difficult to show 

n 

sgn(a) ]J sgn^a) = -1 , (4.9) 

for any a and n. The summations over ji, j 2 ■ ■ ■ and j n can then be performed for 
all a at once, leaving with this minus sign. 

Now we turn to the integrations over Zi (i — 1 ~ n). The convergence on 
the geometric series leads to the successively ordered integrations of z\s for each 
a. We assume here that only the poles at Z{ = z* give rise to terms with physical 
significance. This is the case for three-loop correlators. By simply picking up a pole 
of Zi at — tt — for i = 1 ~ n and using 

1 pi-A(zi;Ai) to 

dzi s ( d \ ( 1 \ 



2th v ' \dki) \ P i-[A](zi-Ai)J 

p dzif(---Zi,---) , £ = 0,1,--- , (4.10) 



d ( d s ' 



d(ad) \dA, 
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we find that eq. fl4.5|) is written as 



where 



R {n) = E Ua7 /•"/ E A*- 2 ([ii - *(ii)D II AKb' -<?')]) 

il \ 41 / ' cr65„ i(^ii) 

+ e E(^n^)^7-7^ ([<i -^ 



H+ m 2 (11,12) *1 

2 n 



l *2 

=n — 2 

x II Oo(b'-^)l) 



™, /■ ■ ■ ^ L/^Ji L/^Ii L/^Ii ^ J 



l+ m 2+ m 3 (11,12,13) *1 



«1 



j'(/ii,*2) 

+ ■ ■ • 

+ E 



(ii |i2 ,.,i„_2) dAj^dAj ^dA ln _ 2 

x H D l ([i j - aiij)]) J] A,([j-^(j)]) • (4-12) 

The integrals in the equation above are with respect to z?'s and we adopt a 
notation 

\i) = z* , . (4.13) 

/ \ n—2 

This expression is in one to one correspondence with the expansion of (Z)i=i x i) 
The number of terms appearing is equal to the number of partitions of (n — 2) objects 
into parts. 

In order to put eqs. ( |4.11|) , (|4.12[) in a simpler form, let us introduce 



mi, m 2 , m n \ 

ii, Hy in ) n 

= - E ^mi([«i - cr(h)])D m2 ([i 2 - <t(i 2 )]) • • -D mn ([i k - a{i n )}) 

<TS<Sn 

= - V l - \ \ . (4.14) 

^ n [<i - [i 2 - <t(z 2 )]^ +1 [i n - a(*„r»+i 1 J 
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In particular 
,-2 

V k, h, 

-3, 1, 0, 
k, k, h, •-, 



n-2, 0, 



= - e A^azi - (7(101) n ^o(b'-^o')]) 

= - £ D ft _ 3 ([ii - t7(ii)])L>!([i 2 - o-(i a )]) 

x II A)(b' -<?')]) 

i(^«i.*2) 

e.t.c (4.15) 



In the [Appendix B| , we prove that 

= if ^m^<n-3 , (4.16) 



mi, m 2) m n 



as well as 

/ mi, m 2 , m fc , 0, 

V k, k, ik, ik+i, 

1 ( mi, • • •, mi - 1, • • • , m k , 0, 



n 

k 

E -j \ - 



n-1 



if J]m £ = n-2 . (4.17) 
In particular, 

f n-2, 0, • ■ ■ \ 1 / n - 3, 0, • ■ ■ \ 1 



Zl ' "' " ' / n [*i - in] 2 \ Zl > 



and 

/ n-3, 1, •■■\ = 1 / n-4, 1, 



1 n b'i-i] 2 

(4.18) 



n-1 



+ ^^f" -3 ' °' "') • (4-19) 

Let us introduce graphs in which the factor — j] 2 is represented by a double 
line linking circle i and circle j to handle the quantities defined by eq. ( f4.14| ) more 
easily. For example, for n=3, 

f 1 °) =^ !_= 2 1 3 (4 20) 
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Using the recursion relation eq. ( |4.17 ) and eq. ( J4.16 ) we have, for n=4, 



2 _ 1 / 2 - 1 
1 2 3 4 ) ~ [zt-ztf \ 1 2 3 



* = ,Jl s ■ (4 ' 21) 



1 1 \ 1 / 1 — 1 1 \ 1 / 1 1-1 



i 2 3 4 ; 4 [ Z * - z * 4 ]2 y i 2 3 ; 3 [** - z 4 *p ^ i 2 3 ; 3 



x 123 + ^ V x 213 
00=0 4 O 00=0 



4 1 2 3 + 4 2 1 3 , (4-22) 



and, for n=5, 



3 \ 1 / 3-1 



1 2 3 4 5 / [zf- z* b f \ 1 234 



2 1 \ 1 / 2- 1 1 




1 2 3 4 5 ; 5 [z\-ztY\ 1 2 3 4 / 4 
1 / 2 1 - 1 \ 



[z* 2 -ztf V 1 2 34 



1 



4 1 2 3 + 4 2 1 3 

(KKK) (XXX) 
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+ 58 2cr\>3 



(4.24) 

From the examples above, it is clear that the graphs for the general case can be 
written down quite easily. 



In terms of the quantities defined by eq. (|4.14 ), we obtain a formula for the 
n-point resolvent : 

where 



W ^ A™" 2 /" f ( n-2, 0, 



'1 



\ - \ - / d \mw g xm 2 /■ /" / mi, m 2 , 0, 
^ . .^AflA, / VflA, / /"' / I ii. i,. 



+ ^ ^ V^A~J V^a-J l^A-J 



=n — 2 

/ mi, m 2 , m 3 , 0, • - -, \ 
V h, k, h, U, • • •, «n / n 

^(i)-(afc)/-/(i &lh<*> 

Here > 1 and the summation (zi, • • • , ik) denotes a set of k unequal integers from 
1, 2, • • • , n and ik+i, • • • , «n in the array represents the remaining integers. 

From eq. ( 4.25Q , eq. ( |4.26| ), we can derive the rule which states how to construct 



the n-resolvent correlator graphically. This rule was shown in sect. 
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4.2 The selection rule in summation and crossing symme- 
try 



In the previous section, we have derive the formula for the n-resolvent correlators, 
which were functions of Q and ji. In order to obtain the loop correlators, we have 
to carry out the inverse Laplace transformation with respect to Q further. In the 
processes, the following functions will play a fundamental role. 

n ( d 
P n (0i,9 2 , ■ ■ ■ 9 n ) = ^2 




11 Bz* 



where 



z* = exp(2r]i cosh ^) , (4.28) 
Pi - = ad = 2r]™ cosh = aM, cosh , ^ = (aM;/2) 1/m (4.29) 
Ai-A* = -(m + lH 2m = -a 2 //; = -(m + l)(aM,/2) 2 , (4.30) 



and 

<9A 



1 cto/ coshm^j cte* / di] \ sinh(m — l)9i 



rj dA sinh mOi <9A j V <9A / sinh 



(4.31) 



The important point here is that the above functions can be written as a sum of 
the product of n factors each of which is a function only of the corresponding 9i and 
/i, so that we can obtain the inverse-Laplace images of P n (9i, • • • , n ) immediately. 
We conjecture that n-resolvent correlator can be expressed in terms of P n > (n' < n). 
In fact, as we will show later explicitly, it is true for the four- and the five- resolvent 
correlator. 

A key manipulation we will use is the partial fraction 

' +T- ^ TT • ( 4 - 32 ) 



[i - - k] [i- k\[k - j] [i - - k] 

One can associate a line from % to j with uzjt- The following identity is responsible 
for expressing P n as a sum of the product of n factors each of which depends only 
on Of 

^ _ 1 / sinh(m — k)a sinh(m — k)(3 \ 

cosh a — cosh (3 \ sinh ma sinh m(3 J 

^rn-k k s j n h( m - j - i + i)o; sinh(m - j - k + i)P 
~^ sinh ma sinh m(3 

(4.33) 
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Let us first work out the cases n = 2, 3, 4 to get a feeling. For n = 2, 

P,M) = ^f = 2 (U) 

2 f ^ f zl^ 'v' V sinh ( m ~ ji ~ ji + sinh(m - ji + h - k x )6 2 
\ dA ) \ Tj J sinh m9\ sinh m# 2 

2 ( Hi) ( —\ y sinh ( m -Ji) g i sinh(m- ji)g 2 , A M , 

\dAj \ i] ) sinhm^! sinhm# 2 



where i\ = k\ — 1. For n = 3, we use eq. Q4.32 ) for the term containing -J?- to 



dA 2 



create a link [1 — 3], which is originally absent. This relates P3 to P 2 . We find 

Pz{0i,e 2 ,6z) 

d V \ h (62,61,771) -h(e 2 ,e 3 ;m) 



2 



<9A/ [1 - 3] 

2 / m— k± ki=l\ /m—k2 k 2 

\dA I \ n ) i , , / \ , , 
sinh(m — j 2 — i 2 + 1)6\ sinh(m — j\ + h — ki)9 2 sinh(m — j 2 + i 2 — ^2)6*3 



A a 4\{4\ I E Ell E El 



sinh md\ sinh m# 2 sinh md 3 

Here fc 2 = Ji +h - 1 = ji- 

This can be repeated for arbitrary n. In the case n = 4, we use the partial 
fraction for the two terms containing and ^jL to create a link [1 — 4], which is 
originally absent. This enables us to relate the case n = 4 to the case n = 3. In 
general, P n is related to P n _i by using the partial fraction for the terms containing 
tj^ 2 - ~ » "~ 1 to create a link [1 — n]. We obtain 

/jt 1 sinh(m - > + - sinh(m - j n _i - a n _i + l)6>i / 4 gg\ 



= 1 



sinh m9£>+i I sinh m#i 



where = je-i + — 1, for £ = 2, 3, ■ • ■ , (ra — 1). Eq. (|4.36|) expresses the 

sih1i('t?7' J\i\0 ' 

Pn(6u ®2, • • • #n) as a sum °f t ne ^-products of the factor — -. Owing to 

sinh mOi 

this property, one can perform the inverse Laplace transform immediately. 

Let us now discuss the restrictions on the summations of 2n — 3 integers ji, i 2 , j 2 , 
■ ■ ■ i n -x,j n -i in eq. (|4.36|) . We write these as a set: 

^n(jl,k,32 ■ ■■in-l,jn-l) 
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{(Ju i 2,j2, ■ ■ -in-ijn-i) I 1 < k < h, 1 < je < m ~ h, for £ = 1, 2, ■ ■ ■ n - 1} 



n-l 



= F 2 {h = 1, ji; fci = 1) IJ nF 2 (i t J e ; k e ) , (4.37) 
where 

Ti{ii,ii\ k e ) = {(ie,je) | 1 < H < ke, 1 < je < m - kg, with k e fixed} . (4.38) 

We will show that these restrictions on the sums are in fact in one-to-one correspon- 
dence with the fusion rules of the unitary minimal models for the diagonal primaries. 
Let us begin with the case n = 3. Define 

Pi = jl +h-H , p 2 = 32 + k 2 ~ 12 , ?3 = h + 12 ~ 1 , , 4 3Q n 

a 12 =pi-l, a 23 =p 2 -l, a 3 i = <?3 - 1 , 

The inequalities on « 2 , J2 are found to be equivalent to the following four inequalities: 

012 + 023-^31 = 2(k 2 -i 2 )>0 . 

012 - «23 + «31 = 2(« 2 - 1) > . 

-ai2 + o 23 + a 3 i = 2(j 2 - 1) > 
au + a 23 + a 3 i = 2(j 2 + k 2 - 2) < 2(m - 2) . (4.40) 



From the third and the fourth equation of eq. ( 4.40 ), the inequality a\ 2 < m — 2 
follows, which is a condition for T 2 ii\ = 1, ji; k\ = 1). Defining a set 

3 

P 3 (a 1 ,a2,a 3 ) = {(ai,a 2 ,a 3 )| ^ a i: - aj > for i = 1 ~ 3 , 

3 

a* = even < 2(m - 2)} , (4.41) 

i=l 

we state eq. ( [4.40|) as 

•^(ji,^, J2) = ^3(012,023,031) • (4.42) 

We also write 

^20i)=^2(n = l,ii;A; 1 = l) = P 2 (a 1 2) . (4.43) 
for the case n = 2. 

Eq. Q4.41|) is nothing but the condition that a triangle be formed which is made 
out of ai, a.2 and a 3 and whose circumference is less than or equal to 2(m — 2). It is 
also the selection rule for the three point function of the diagonal primaries in m-th 
minimal unitary conformal field theory |i~7 |. 
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For the case n = 4, introduce p^ = jz+k$—iz , 034 = p 3 — 1 , g 4 = ,73+23 — 1 , 041 = 
g 4 — 1 . We find 

^2(^3, h\ h) =^3(031,034,041) ( 4 - 44 ) 

The restrictions on the sum in the case n = 4 can be understood as gluing the two 
triangles: 

^4(ji,*2, h) = ^3(012, 023, o 3 i) nP 3 (a 34 ,a 4 i,a 3 i) 

= ^4(012,023,034, 041; o 3 i) . (4.45) 

The allowed integers on a 3 i are naturally interpreted as permissible quantum num- 
bers flowing through an intermediate channel. As one can imagine, eq. ( |4.45| ) is not 
the only way to represent the restriction: one can also represent it as 

^4(ji,«2, J2,«3, J3) = ^3(012,024,041) n£>3(a23,a 3 4,a 2 4) 
= ©4(012, ct23, 034,041 ;a 24 ) , (4.46) 

which embodies the crossing symmetric property of the amplitude. 

The restrictions in the general case n are understood as attaching a triangle to 
the case (n — 1). To see this, define 

(4.47) 



Pi = 


je + kt-it , qe = 


+ - 


1 


' for £=l,2,---n 

1 


ae,i ~- 


= Qe - 1 , 


0^,£+l = Pi — 


1 


1 < i n 


-1 < fcn-l, 1 < jn 


_i < m — 


1, 


we derive 




On— l,n + Q"n,l O n _ 


-1,1 = 2(j n _ 


-1 


-i)>o, 




O-n— l,n O n i + 


-1,1 = 2(fc n . 


-1 


- i n -l) > , 




~Q"n—l,n + O n i + fl n ,- 


-1,1 = 2(z n _ 


1 


-i)>o, 




On-l,n + O n) l + a n _ 


-1,1 = 2(j n _ 


-1 


+ A;„_i - 2) < 2(m - 2; 



(4.48) 

The restriction on i n _i and j n _i are, therefore, D 3 (a n _i ira , a W) i, On-i^), which is what 
we wanted to see. All in all, we find 

, ' ' ' l n-liJn-l) 

= ^(On-l.n, O nj i, O n -l,l) H T n -l(jl, ili J2i ' ' ' *n-2 3 Jn-2) 

= ^3(On-l,n, O n ,l, Qn-l,l) H XV-l(°l,2, 02,3, ' ' ' O n _ 2 ,n-l, On-l,l5 a 3,l, 04,1, • • • a n -2,l) 

= P n (oi,2, 02,3, ■•• o n _i jn , a n> i; 03,1, 04,1, "■■ £*n-l,l) (4.49) 

From now on, a shortened notation P n (ai j2 , a 2 , 3 , • • • a n -i,n, o n ,i) is understood to 
represent P n (oi,2, o 2 ,3, • ■ • a n -i,n, o„,i; a 3j i, a 4i i, ■ • • a n _i,i). 
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Putting eq ( |4.36| ) and eq. ( |4.49| ) together, we obtain a formula 



' drf 



-1 



^ En 



sinh(m — kj — l)9j \ sinh(m — k\ — 1)9\ 



T> n \i=2 



sinh mOj 



sinh mOi 



(4.50) 



where T> n means T) n {k\ — 1, • • • , k n — 1). Once again, the fact that the different 
divisions of T> n into n — 2 triangles are embodied by this single expression is precisely 
the statement of the old duality. 

The object P n (6 l 1 , 9 2 , ■ ■ ■ 9 n ) is equipped with 9j and kj for j = 1, 2, ■ • ■ , n and any 
V 3 (ki — 1, k 2 — 1, k 3 — 1) obeys the rule of the triangle specified above. It is, therefore, 
natural to visualize this as a vertex which connects n external legs corresponding to 
n loops. The vertex can be regarded as a dual graph of an n-gon that corresponds 
V n (kx - 1, • • •, k n - 1). 

Due to the formula ( |3.30| ), we can obtain the inverse-Laplace image of eq. ( |L50| ) 
with respect to Q (i = 1, ~ n) immediately: 



£-Mp„ w ,.,,„)]=-(^)(^)"(f)"E[n^ & (^ 



1 /flMv 



m(m + 1) 



. (4.51) 



4.3 Four-loop correlators 

In this subsection and in the next one, we show how to perform the inverse Laplace 
transformation of the resolvent correlators to get loop correlators in terms of loop 
lengths in the case of n = 4, 5. It is necessary to put 

R^(9 1 ,...,e n )^R^(z* 1 ,...,z:)\ At=A (4.52) 

in a manageable form to the inverse Laplace transform. Let us recall that P n (6i, ■ ■ ■ , 9 n ) 
can be inverse Laplace transformed immediately. If R^ n \9i, • • • , n ) is expressed as 
a polynomial of Pj{9\, • • • , 9 n ) and their derivatives with respect to A, the inverse 
Laplace transform can be carried out immediately. In fact, this turns out be true 
for n = 4 and 5, which we will show explicitly in the following. 

The important point we will use in the following is the fact that when one of the 
loops shrinks and the loop length goes to zero, the n-loop amplitude must become 
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proportional to the derivative of the (n — l)-loop amplitude with respect to the 
cosmological constant. We represent this fact by 







(w+ih) ■ ■ ■ w + {Q) -> oc — {w + {i 1 ) ■ ■ ■ 



(4.53) 



when n-th loop shrinks. Then in the limit, PS n '(9i, • ■ ■ , 9 n ) must satisfy the following 
relation: 

'R {n) (9 l ,---,9 r ;' ' ' ° 



C 



-l 



oc C~ l 



—R^(9 1 ,---,9 n - 1 ) 



(4.54) 



This relation restricts the possible form of R^ n \9i, ■ ■ ■ ,9 n ) . We should note here 
that the inverse Laplace image of P n (0i, ■ ■ ■ , 9 n ) satisfy the following relation 



£ 1 [Pn(9l, • ■ ■ ,9 nJ 



This fact follows from substituting 



K 1 _k A (M£ Tt 



'ML 



Mt n < 1 



(4.55) 



(4.56) 



in eq. (|4.51| ) and picking up only the k n = 1 parts. 

Because we want R^ n \9i, ■ ■ ■ , 9 n ) to be expressed as a polynomial of Pj and their 
derivatives with respect to A, we need here to introduce a notation 



[^Pl,-M^ > 32,-,h+ii-l ' ' ' Pn-ig+l,-,n) ($1, $2, " " ' 9 n ) = 
1 



(4.57) 



m 



^ Pil(9cr(l), • • ■ , 9 a (i 1 ))Pi 2 (9a(j 2 ), ■ ■ ■ , ^(ja+ia-l)) ■ ■ ■ Pi e {9 a ( n -i e+ i), ■ ■ ■ , CT ( n )) , 



<reVn 



where V n represents the permutations of (1, 2, ■ • • , n). To be more specific, for 
example 



1 



[<S-Pl23-P234] (01, #2, #3, 6> 4 ) = — ^2 P3(9 a (i),9 a!y 2),9 CT ^)P 3 (9 a (^2),9a(3),9 a 



(4), 



0-GP4 



[«5Pl234-P34] (01, 9 2 , 9 3 , 64) = — P±{6(t(1)i 0(7(2), 0(7(3) 0ct(4) )-P2(0<j(3), 0(7(4)) ■ (4.58) 



It is convenient to represent P n (9i, ■ ■ ■ ,9 n ) by an n- vertex which connects n 
external legs. For example for n = 2, 3 and 4 



P 2 (0i,0 2 ) 




^3(01,02,03 
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Pa{0i, 02,03, 9 a) 




(4.59) 



The n- vertex can be regarded as a dual graph of the n-gon (polygon) which corre- 
sponds to V n . In terms of these vertices, let us express eq. (|4.58|) as follows. 



[5P123-P234] 02, 03, 0a) 





[SPi23aP3a] (01,02,03,0A) 




(4.60) 



The relation eq. ( 4.55| ) can be represented, for example, as 




OC 




(4.61) 



in the case of n=5. 

Now we are concerned with the case of n=4 first. Let us recall that for n=3 

R {3) (01,02,03) ocP 3 (0i, 02,03) ■ (4.62) 

(0i, 02, 03, 0a) must include a term which becomes proportional to ^3(^1,^2,^3) 
in the limit M£ 4 — > 0, which is P±(0\, 02, 03, 0a)- R^'(0i, 02, 03, 0a) ma y also include 
terms which vanish in this limit. Such terms must consist of the product of two 
multi-vertices which have 6 external legs in total. 

By explicit computation, we find 

^R {4) (01, 02,03,0A) |= ^Pa(01, 02, 03,0A) ~ [SP123P23A} (01,02,03,0A) 
+ [SP 12U P U ](0 1 ,02,0 3 ,0 4 ) 
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where the prime represents the differentiation with respect to A. 

By performing the inverse-Laplace transformation and renormalizing, we obtain 
the complete answer for the macroscopic four loop correlator: i 



V(4 



W +(Q) = A{ uslon (£ u ■■■£,) + AT ldual (£ u ■■■£,) , 



(4.64) 



where 



and 



Here 



mlm 



1 4 d 



M X 2(l-i) 

2 



eii^ !±( m i 

T> 4 j=l 



(4.65) 



^residual / £ 



4! 



E 



l) 



m(m + 1) 



i=i 



eeii [ 5 . 

x> 3 i=i 



123 5 /234 



+EEn[^. 

©4 3=1 



1234 ~ R /34 



M£ aU) 



(4.66) 



1234 
5f 



5 



(-K'l-fex/m) ^1 

/m ) j 

{Kl-ki/m, Kl —kz/mi —h/±/r, 

k'./ra) • 



(4.67) 
(4.68) 
(4.69) 
(4.70) 



3 This was briefly reported in 
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We have introduced V 3 = V 3 (kt -l,k 2 - 1, k 3 - 1), V 3 = V 3 (k[ -l,k' 2 - 1, k' 3 - 1), 
D 4 = T>4(ki — l,k 2 — 1, &3 — 1, A; 4 — 1) and X>g = T> 2 (k[ — l,k' 2 — 1) and have defined 
the convolution [A © B] (Mi) with respect to £ by 

[A © 5] (Mi') = [ d£'A(M£')B(M(£ - £')) , (4.71) 
Jo 

in particular 

[A 5] (M£) = f d£A(M£')5(£ - £') = A(M£) . (4.72) 
Jo 

The important point to note is that J^^ stdual seems to represent the contribution 
from loops with mixed momenta, that is, a single loop seems to have two distinct 
parts each of which have different momentum. If two of the loops l\ and £ 2 touch 
each other on two point, the two-dimensional surface break into two surfaces and the 
loops £\ and £ 2 also split into two pieces respectively. We infer that this configuration 
of two-dimensional surface may have connection with A? 4 esidual . For m = 2 (pure 
gravity), J^ sldual vanishes. This fact is consistent with the above consideration 
because in pure gravity there is no loop with mixed boundary condition. 



4.4 Five-loop correlators 

Let us now turn to the n=5 case. R ( - 5 \9 1 , ■ ■ ■ , 9 5 ) include a term which is proportional 
to 




(4.73) 



corresponding to the first term in eq. (|4.63|) . Corresponding to the second term in 
eq. (|4.63j ), R^(9i, ■ ■ ■ , 9 5 ) must include a term which consists of the product of two 
multi-vertices with 7 external legs in total. The possible form is 



d 

<S7^-(Pl2345)-P45 



r 9 

+ b <SPi2345T^-(-P45, 



+C 



d 



r 9 

+ d S ^1234 7^(-P345) 
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(4.74) 



In the limit Mi* — ► 0, it becomes 







(3a + c) [S—(P 12U )P M \ + (3b + d) 







+2c 



.9 



2d 



5Pl234 7^-(-P34y 
d 



(4.75) 



We require this expression to be proportional to the A-derivative of the second term 



in eq . (g6g ). We find 



(l-c)/3, 6=(2 + c)/3, d 



-1 - c 



(4.76) 



R^(9i, ■ ■ ■ ,65) may include terms which vanish in the limit under consideration as 
well. They must consist of the products of three multi-vertices with 9 external legs 
in total. As one of the such terms we have 



[1SP1234-P34-P45] — 2[5Pi234-P23-P345] + [5-Pl23-Pl24-P 



235 J 






(4.77) 



( There are some other combinations which satisfy the conditions. ) R^ 5 '(9i, • • • , #5) 
must be expressed as a linear combination of the above three types of terms eq. ( |4.73j ), 
eq. ( |4.74|) and eq. ([4.77 ) if the assumption under consideration is true. By explicit 
calculation, we have found in fact that R^(9i, ■ ■ ■ , 9 5 ) can be expressed as a linear 
combination of eq. Q4.73D , eq. ( |4.74j) and eq. Q4.77Q : 



R( 5 '(9i, 9 2 ,9 Z) 9±, 9 5 



dA 



+ 



SR 



12345 



2 dA +3 dA 



P${9i } 9 2} 9s, 64, 9 5 ) 
P45 (#i, #2, #3, #4, #5) 
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SP,2-,u l Sr + 4:Sr I P 



:>> 15 



h, #2, &3, #4, #5 



OA dA 

+ [(2P1234P34P45 - 4P 123 P2 3J P345 + 2P 123 P 12i P 235 )] (9 1 , 9 2 , 9 3 , 9 4 , 6 5 ) . (4.78) 

Following the same procedure as obtaining eq. ( (4.64 ), we find the complete an- 
swer for the five loop amplitude: 

V(4)---™ + (4)! 



A fusion^ £i . . . 4) + ^rwMual-l^ . . . ^ 

Al esidual - 2 {ix---h) , 



(4.79) 



where 



m(m + 1) jj[ ^dfi 

residual— 1 
1 



Z © 5 J=l 



(4.80) 



^ esidMi_i (^i---4) 



m(m + 1) 



5! 



E 



X 



EEn[^ i2345 pf] (M4 -)) 



£>S j=l 

EEn[^ 1234 ©^ 345 ] wo-)) 

2?4 I>' q i=l 



(4.81) 



ik UtE v 



Mx3(l-i) 



•2 E E E f[ [£f 345 © © *T] (M£ a(j) ) 
v 5 v> 2 t>2 j=i 

+4 E E E f[ [£f 34 © Bf Pf 45 ] (M£ a{j) ) 

T>4 V' 2 V'-l j=l 

- 2 E E E n [ fil23 © 5 f 4 ® s f 35 ] W*>) 



(4.82) 



T>3 V' V" 3=1 
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Here (J^)^^ means the derivative acting only on the left(right) part of the con- 
volutions. The rest of the notations here are similar to those of the n = 4 case and 
will be self-explanatory. 

We conjecture that can be represented as a polynomial of Pj, j < n and their 
[i derivatives: the final answer would then be obtained by convolutions of various 
Bj's and their derivatives. If the conjecture is true in fact, the power counting 
argument tells us that the j-th term of turns out to be represented by a figure 
which consists of the products of j multi- vertices with n + 2j external legs in total. 
We hope that, for higher loops, R^ n \9%, • ■ ■ , 9 n ) can be put in principle in a form 
such as eqs. ( |4.64j ), fl4.79|) in the same manner as we have determined the four and 
five loops from the lower ones. From eqs. Q4.64 ) ~ ( |4.65 ) and eqs. ( 4.79|) ~ (|4.80| ), 
we guess that the multi-loop correlators, in general, would include the following 
term corresponding P n (9i, • • • , 9 n ) : 



= !— fr<f-) 

m{m + 1) M 3 ^d\i> 



8 \n-3 



M\-M)(i-^) 
2 



(4.83) 



4.5 Four-point functions from loop correlators 

The four-loop correlators we found in eqs. (4.64) ~ ( |4.66| ) do not diverge when the 
loop lengths approach zero, so that we expect that the loop operator can be replaced 
by the sum of the local operators in this case. Let us derive the four-point functions 
of the scaling operators, applying the expansion of loop to the four loop correlators. 

First, consider the part A{ us lon (£i, ^2,^3,^4) and expand this in terms of the 
modified Bessel functions I v {ti). From eq. (|4.65|) we have 



1 



m 2 (m + l)\2j 2 31 

- 2 - f^V 2 " -LJLIy,~[1—K i (Mi,) } . (4.84) 
m(m + l)H2 j MdMy^J}, 2 W l 3> \ K ' 

Since we can prove the relations 

-K^ p (z) = J2 (P + 2n ) I]p+2n\(z) , (4.85) 
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and 



E (P + 2n) (p + 2n(p + n)) I\ p+ 2n\ (z) 



for < p < 1, we obtain the following expression for A{ uswn : 



(4.86) 



m(m + l) 2 V 2 / 



*eii e 

X>4 j = l \%*=— OOy 
4 



1 1 4 

3=1 



m 



m m 



(4.87) 



Comparing eq. ( |4.87| ) to the expansion of loop operators ( |3.65| ), we obtain the 



following contribution of the four-point functions from A{ uslon (£i 



fusion 



> • " " j 

4 , fci 



\j=i 



a 



k\k2kzk± ~< 



2m 2 (m + l) 2 V 2 / 

4 

2(m+l)-E[% + 2 ^' + 



where 



5=1 

m— 1 

Ck 1 k 2 k 3 k4 = E/ Ck^k'Ck'kski 
k'=l 



(4.88) 



(4.89) 



Next, let us consider >/ 44 esidua '(£ 1 , ■ ■ ■ ,£4) part. We can prove that the convolution 
of two modified Bessel functions K p (M£) and K p i(M£) is expanded in terms of 
I U (M£) as 



ifp K p , 



{Mi) 



' — ' E {«(2n + p + p')/|2nWlW 

n=— 00 



- n{2n + p - p') I\ 2n+P - P >\ (M£) } 



(4.90) 



The above relation is easily derived in the space of Laplace transformed coordinate 
using the relations 



C [ I P {M£)\ 



^ e -pm9 

M sinh m9 



for p > — 1 



(4.91) 
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and 



K P (M£) 



C [ I- P {M£) - I P {M£)\ 
2 sinhpm^ 

for p > — 1 



M sinhm^ 

because the convolution in £-space corresponds to the product in ^-space. 



(4.92) 



From eqs. ( 4.85| ) and ( |4.9C ) we obtain the following expression for the residual 
part of the four loop correlator: 



A 



residual 







(M 



1 

4! 



E 




m 2 {m + l) 2 \ 2 

EE(-^ 

lk 2 k 3 C jfc^fc 4 + Ck 1 k 2 k3k 4 ,Ck' 1 k' 2 ) 



£(±)ni(2m + (l- 



x£(±)n 2 (2n 2 

+- 

x(2n 3 + 1 ' 



±fl 



m 
m 



|(2n 2 +(l-£)±(l- 



,(M4 



(2), 



1(2-3 + ^)1 



M (M4( 3) ) (2n 4 + 



m 



l(2«4+^) 



*4\|(-W^<x(4), 



(4.93) 



where C, 



kk< 



>kk' 



From the above expression and the expansion of the loop operator ( p. 65 ), we can 
obtain the contribution of the four-point functions from j\^ esidual part. The explicit 



expression, however, would be complicated. 

Let us comment on the role of the coefficients 



k2k^k\k' 2 k A 



in eq. ( [4.93] ). At first sight, it appears that we would have l v [Mi^) with integer order 
in eq. ( [4.931 ). For example, in the case of k\ = k[ or k\ + k[ = m, we have l v [Ml\) 
with integer order. Terms including l v [Ml^) with integer order as a factor cannot 
be explained from the viewpoint of the local operators. These terms, however, are 
cancelled due to the coefficients {—C^^Cy^y^ +Ck 1 k 2 k- i k i Ck' 1 k' 2 )i so that we do not 
have lyiMii) with integer order in eq. ([4.93j) after all. 
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5 Summary and discussion 



In this paper we have investigated the correlators of macroscopic loops and those 
of local operators in the unitary minimal models coupled to two-dimensional gravity 
using the two-matrix model. We calculated the general multi-resolvent correlators, 
and examined one- to five-loop correlators explicitly. 

From these loop correlators we obtained the correlator of the scaling operators 
by applying the idea |3(| that the macroscopic loop can be replaced by a sum of 
local operators, to the case of the two-matrix model. We found that there exist the 
fusion rules for the three-loop correlators, which are similar to those for the three- 
point functions of the gravitational primaries. From the three-loop correlators, we 
deduced the three-point functions of the scaling operators, and found the gravita- 
tional descendants as well as the gravitational primaries satisfy the fusion rules of 
the same kind. These fusion rules for the loops can be considered to express those 
for all of the scaling operators in a compact form. 

At the (m + l,m) critical point in two-matrix models, the scaling operators 
dj, j = (mod m + 1) have no counterparts in the BRST cohomology of Liouville 
theory coupled to the corresponding conformal matter. In [33], these operators were 
argued to be boundary operators which couple to loops in the case of the one-matrix 
model. It was also shown explicitly that one of them, corresponding to d m +i in the 
case of the unitary matter, is a operator which measures the total length of the 
loops. 

We examined the role of the rest of these operators. We showed, in some exam- 
ples, that the operator B n couples to the points to which n parts of several loops 
are stuck each other. In other words, the operator B n connects n parts of loops 
together. We think these operators play an important role concerning the touching 
of the macroscopic loops. The emergence of these operators in matrix models can 
easily be understood from the viewpoint of macroscopic loops and their expansion 
in terms of local operators. 

In sect. (| we examined the property of the multi-loop correlators. We pointed 
out that the structure similar to those of the crossing symmetry in the underlying 
conformal field theory can be seen in the loop correlators. This structure appears 
in the selection rules for the summations in the expression of loop correlators. 

We calculated explicitly four- and five-loop correlators. From the expression of 
these correlators, we inferred that these include the contribution from the loops with 
boundary condition specified by more than one momentum. We guess this property 
can be understood as follows. When two loops touch each other on two points, each 
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loop breaks into two pieces. Since a single loop breaks into two pieces, the broken 
pieces can have distinct momenta. The configuration probably have non-vanishing 
contribution to the amplitude in the case of the four- and five-loop amplitudes. Note 
that matter degrees of freedom are fixed only inside the loops when we calculate the 
loop correlators. Each loop therefore represents a superposition of loops with various 
momenta. Distinct lattice elements on a single loop can have distinct momenta. 
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Appendix A 



Here we collect some formulae concerning the modified Bessel functions. I v {z) 
and K v (z) are linearly independent solutions of the Bessel equation 



( 9 \ 2 2 2 



Z v (z) = . 



Ifj,(z) can be expanded as 



■Z\2n 

L) 



and K„(z) is defined as 



vr/2 



[/_„(*) - /„(*)] . 



Sin Z/7T 

We collect another useful formulae in the following. 

- — K u (z) = K v _ x (z) - K v+l {z) 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(^)V^)] = (-l) n ^-^-n(^) 



(A.5) 



TO) .i r ^ (f)exp( _i_£!±c!) 

[Re z > 0, i?e C > 0, |arg(z + C)| < tt/4] (A.6) 
Introducing the Laplace transformation of a function f(£) of loop length £ by 



/(*) = / dl e-Vf(x) , 
J Jo 



(A.7) 



we have the following relations on the Laplace transformations of the Bessel func- 
tions: 

I e -pm8 



I P (M£) 
K P (M£) 



M sinh mO 
2 sinh \p\m9 



P> -1, 



M sinh m# 

~ 1 \p\K p (M£)\ = 2coshpm6, 



where ( is parametrized as 



C = M cosh 7716* , 



and we introduced a notation 



K P {M£) = ^}^K p (M£) . 

7T/ Z 



(A.8) 

(A.9) 
(A.10) 

(A.11) 
(A.12) 
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Appendix B 



In this appendix, we prove the recursion relations for 

with ^2 m £ < n — 2 



m x , m 2 , m 3 , 

k, k, «3, 



(B.l) 



introduced in the text. The proof goes by mathematical inductions. 
We will first prove the simplest case 



m, 0, 
k, k 







1 



for m = n — 2 



fB.2l 



n n [<i 



and 



m, 0, 
k, k 



for m < n — 3 



(B.3) 



Assume that eq. ( p.2| ) and eq. (|B.3|) are true at n. Without loss of generality, 
%x can be taken to be 1. Let us consider the left hand side of eq. ( |B.2|) or eq. ( p.3| ) 
in which n is replaced by n + 1. To compute them we observe that the elements 
of S n+ \ are generated by associating n different ways of inserting [n + 1] with each 
element a G «S n . In the case where [n + 1] is inserted in between [1] and [c(l)], this 
contribution is equal to 

1 [1 — o-(l)]" 1 n 1 

~ t 1 - t 1 - (« + + 1) - *(i)] Jk [j - ' (R4) 

The contributions from the sum of the remaining n — 1 insertions are found to be 
equal to 



1 



_5 \l-a(l 



1 



„es n r Al)] m HI) ~ (n+ 1)][1 - (n+ 1)] [j - 
Here we have used 

1 1/1 1 



1 



[j — (n + 1)] [{n + 1) — m] [j — m] \ [j — (n + 1)] [m — (n + 1)] 



n-l 



Note also that - = IlVb' ~ cr(j)]- Putting eqs. 

J'=l J'=2 

flB.5|) together, we find 



(B.5) 



(B.6) 



and 



m, 
1. 



0. 

2. 





n + 1 



n+l 
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t 1 " ^(!)] m K 1 ) - (» + W - (» + !)] M b " *0")] 
Factorizing the expression inside the bracket { ■ ■ • }, we have 

/ m, 0, • ■ ■ , \ ^ f m-l, 0, \ 1 

(B.8) 

Then from the assumption, eq. ( |B.2| ) and eq. ( |B.3| ) are also satisfied when n is 
replaced by n + 1. On the other hand for n = 3 eq. ( |B.2| ) and eq. ( |B.3| ) are clearly 
true, so we have proven the relations. 

Now we turn to the more general case the proof of which is a straightforward 
generalization of the one given above. To derive the relations 



mi, m k , 0, 

^1; ' ' ' Ik) ; ^Tl 



k 

, for ^m^<n-3 , (B.9) 



and 



mi, m k , 0, ■-, 

I'll ' ' ' Iki " "j In ) 

_ ( mi, rrij-1, m fc , 0, \ 

j^l V a "l» "' *i» '"' "' /n-lb _n l 2 

for J]m^ = n-2 . (B.10) 



Let us assume eq. ( |B.9| ) at n . 

We take i£ = £, I = 1 ~ without loss of generality. The way in which the 
elements of S n+ % are generated is the same as the one given above. In the case where 
[n + 1] is inserted in between [£] and [a (£)} £ — 1 ~ k, the contribution is 

_ v 1 1 \l - a(£)} me 



creSn 



1 1 ^ 1 



t n t--^t • ( B - n ) 



[(£ + 1) - a(£ + [k - cr(A;)] T "*+ 1 .^...^ [j - a(j)\ 

The contributions from the sum of the remaining insertions are 



n-l 1 j 



i(i i 



X 



ifeiL) ^'^ - + [a^(l)-(n + l)]; ' 

(B.12) 
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Here pe £ = 1 ~ k are such that a Pe (l) = i. Using eq. ( |B.6j ) again, we find that 
this equals 

k 1 1 1 

~ .SnS t 1 - ' " P - <^)] TO < - (n + 1)][^ - (n + 1)] 

1 1 n_1 1 

P + 1) - <r(f + " > - ? J]..,. [; - <x[j]] • (R13) 

Putting eqs. ( |B.11| ) and ( |B.13| ) together, we find 

( m x , m k , 0, • ■ • , ^ 

1, •■• k, k + 1, n + 1 



n+ 1 



\ - \ - 1 1 l 1 V^-(n+l)]; J 

* i*. I 1 - ^(l)] mi+1 [* - ^)) mt W) -n}[£- n] 
1 1 n_1 1 

X P + 1) _ a( £ + l)]m 4+1+ l ■■■[*_ flr(fc )]m fc+ l .JJ^ g - ' (B ' 14) 

Factorizing the expression inside the bracket, we have 

f mi, m k , 0, \ 

I 1, ••• k, k + 1 n+l) n+1 

= A^/mi, rrij-l, m k , 0, \ 1 

(B.15) 

Then from the assumption eq. ( p.9| ) at n, eq. (P-9| ) in which n is replaced by n + 1 
is also satisfied. On the other hand for n = 3 eq. ( P-9|) is clearly satisfied, so we 
have proven eq. (|B.9|) and eq. (|B.10|) . 
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